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“The Canadian Journal of Chemical Engineering” 
formerly “The Canadian Journal of Technology” 


The National Research Council of Canada and The 

Chemical Institute of Canada have entered into an agreement 

| ow hereby the N.R.C. publication “The Canadian Journal of 
Technology” was transferred to the Institute in 1957. 


Commencing with this issue, the journal will be pub- 
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the Institute in Edmonton, Montreal, Saguenay and Van- 
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a lished by the Institute in a new format. The name of the 
journal has been changed to “The Canadian Journal of 
Chemical Engineering”, and it will be adapted to the interests 
of Canadian chemical engineers, carrying both research and 
review papers. 

This new journal will provide a Canadian medium for the publication of progress in chemical 
engineering and allied technical fields in this country, and, with the full support of members of the 
profession, will be a valuable contribution to the advancement of this field of science in our country. 

; Lobfford Purves 
Director of Publicity 
A NEW ENDEAVOUR and Publications 
H. P. Godard 
te conversion and renaming of the long established couver. Chemical engineering student groups are active at 
Canadian Journal of Technology is designed to provide the University of Alberta and the University of Toronto. 
additional service to the chemical engineering members “Chemistry in Canada” has been the news journal of 
= of The Chemical Institute of Canada, and to provide a the Institute and will continue to report the general news 
national publication forum for the increasing amount of of interest to Canadian chemical engineers. The Canadian 
chemical engineering research and development being Journal of Chemcial Engineering will be strictly 
ute done in this county. It recognizes the growth of the technical. Its policy is described on the following page 
ons chemical —— pee in Canada Lage the increas- by editor W. M. Campbell. 
. ing proportion of chemical engineers i S ! ee . a . 
ng. [or . engineers in the Institute. A. Credit for the idea and timing that have led to this 
ent glance at the registration figures for chemists and chemical ed 7 . ; a 
5 ; ees P a new venture goes to past president Roger Gaudry (1955- 
engineers now in our universities makes it clear that this e ; : : : : 
EE ; 1956) who pushed the idea so energetically and enthusi- 
step 1S being taken none too soon. ° « ° 7 - ° 
ind : ; : : astically that in short order the idea has become a reality. 
The Institute has always had a basic policy of serving ihe ch thi wae anita . 
ed, the chemical engineer as well as the chemist. The eg. Wa oe Sot Ch — 7 ae “Chain 
Chemical Engineering Division is one of the largest and the “Canadian Journal of Chemistry” or the ana 
i most active of the eleven subject divisions in the Institute. Journal of Biochemistry and Phy siology” is being made 
: It has held annual divisional meetings since 1951 and par- available to chemists without charge so that at long last 
ticipated strongly in the annual conventions of the Insti- all members receive both a news journal and a technical 
sail tute. It met jointly with the A.I.Ch.E. in Toronto in 1953 journal as part of the services rendered by the Institute. 
a i and is planning another joint meeting in Montreal in 1958. The time is ripe, the immediate outlook for papers is 
as Local chemical engineering groups are active within good, now it depends on you-—the readers. Let’s all help 


whenever possible. Tuum Est! 
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The First Issue 


DITORIALS will not be a regular feature of The 

Canadian Journal of Chemical Engineering but 
will appear only on special occasions. The first 
issue is undoubtedly a special occasion. I would 
like here to review the policies which have been 
laid down for the journal and to ask for help and 
suggestions from the readers, writers and reviewers. 


C.J.Ch.E. will be published every two months by 
The Chemical Institute of Canada. It is a continu- 
ation of the “Canadian Journal of Technology” 
formerly published by the National Research Council, 
Ottawa, and it will maintain the high technical 
standard of this former journal. The essential differ- 
ence, however, is that C.J.Ch.E. will be devoted 
entirely to chemical engineering and closely related 
subjects. For example, applied chemistry is a closely 
related subject if it is slanted towards chemical 
engineering. The dividing line here is somewhat 
indefinite but it should become clearer after the 
first few issues have been published. 

The basic aim of C.J.Ch.E. is to provide a service 
to all chemical engineers in the country. With this 
objective before them the editorial board has recom- 
mended that each issue should have some research 
and some non-research papers with a yearly aver- 
age of about 60% research and 40% non-research. 
The non-research papers will deal with subjects such 
as recent developments in unit operations, flow 
sheets and descriptions of new Canadian processes, 
waste treatment and pollution problems, equipment 
experiences, economics, operations research, and 
technical notes. Papers of general interest to chem- 
ists as well as to chemical engineers will be pub- 
lished by Chemistry In Canada as in the past. 


Publication is essential for the research man, so 
we do not anticipate any problem in getting re- 


t 


W. M. CAMPBELL, Editor. 


search papers, provided the technical standard of 
the journal is high. But publication is less essential 
for the non-research man—the plant engineer, the 
design engineer, etc.—and yet he is interested in 
finding out what his colleagues in other plants are 
doing. In turn he must be prepared to write about 
the things he has learned and the place to do this 
is in C.J.Ch.E. Thus, this is an appeal to all chemical 
engineers, whether engaged in teaching, research, 
plant operations or design, to contribute to the 
journal and, what is just as important, to write and 
tell me what kind of articles we should publish. 
The more specific these suggestions are the more 
help they will be. 

There is another group of men who are essential 
to the success of the journal—the reviewers. We 
expect to publish about 50 papers per year and 
since each paper is reviewed by two men we will 
need about 100 reviews per year. Quite a number 
of men have volunteered to help and over 30 have 
already reviewed papers. These reviews have been 
well done and I would like to thank all these men 
for the time spent on studying the papers and for 
the helpful suggestions which they made. But we 
do not want to overwork these men, so I am 
appealing for more volunteers, particularly those 
who are familiar with recent developments in unit 
operations. Just drop me a line giving the field 
with which you are most familiar. 

We have an excellent start on this new endeavour 
—many good papers have been submitted and 
many more have been promised. The success of 
the journal will depend on three groups, the authors, 
reviewers and readers. It is only with the continued 
help of the first two groups and the expected help 
of the third group that the C.J.Ch.E. will become 
the success that it now promises to be. 
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Effect of Gas Thermal Conductivity 


on Local Heat Transfer in a 


Fluidized Bed’ 


A, JACOB? and G. L. OSBERG* 


Bed-wall heat transfer measurements in fluidized beds, as well as 
experiments with various heater probes, indicate that the smaller the 
heat transferring surface the larger the heat transfer coefficient. Further- 
more, temperature gradients in fluidized beds are very small, except near 
the hot surface. It is, therefore generally believed that particle-particle 
heat transfer coefficients in the bulk of the bed are very large, probably 
because of the transient character of the heat transfer mechanism. To 
shed more light on this mechanism the authors measured heat transfer 
coefficients in fluidized beds of glass beads with a horizontal hot wire 
probe. The diameter of the wire (0.0052 in.) was of the same order of 
magnitude as the particle diameters (0.0115, 0.00602, 0.00242, 0.00122 
in.). Eight gases with thermal conductivities ranging from 0.0048 to 
0.097 btu(ft.) Chr.) (°F.) at 32°F. have been investigated. A heat 
transfer mechanism based on transient heat flow is proposed. The ex- 
perimental data were correlated by the formula 


h = ho (1 — €) (1 — e7 Pkt) 


Where h 
ho,p, are empirical constants 
€ void fraction 


is the observed heat transfer coefficient 


ky the thermal conductivity of the gas 


Introduction 


characteristic feature of a fluid- 
A ized bed in contact with a hot 
surface, e.g. the column wall, is the 
absence of temperature gradients in 
the bulk of the bed. Such gradients 
exist, however, near the hot surface. 
Hence, it is there that most of the 


ticles giving rise to higher heat trans- 
fer coefficients in fluidized beds than 
in still beds or empty tubes. Accord- 
ing to van Heerden and co-workers 
‘) a laminar film of the particle- 
fluid suspension which is formed be- 
tween heater wall and fluidized bed 
constitutes the heat controlling ele- 


resistance to heat transfer occurs. 


To account for the high thermal 
resistance near the hot surface Leva 
and Grummer “!), Dow and Jakob (?) 
and Wicke ‘3) use the concept of a 
laminar gas film at the heater surface 
which controls the heat transfer by 
steady-state conduction. The thick- 
ness of this film is decreased by the 
scouring action of the moving par- 


ment in the case of large heater sur- 
faces. Beyond this film the high eddy 
conductivity of the bed results in a 
sharp decrease of the thermal re- 
sistance. The same authors were able 
to show that the heat transfer co- 
efficients depend on the size of the 
heat transfer surface, high coefficients 
being obtained with small surfaces. 
Experiments with small heater probes 


1 Manuscript received February 2, 1957. 


Contribution from Division of Applied Chemistry, National Research Council, Ottawa, Ont. 
2National Research Council of Canada Postdoctorate Fellow, 1955-56; now with Ciba Limited, 


Switzerland. 


3Head, Chemical Engineering Section, Applied Chemistry Division, N.R.C., Ottawa, Ont. 
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performed by Mickley and Fair- 
banks (5) and Osberg ‘®) confirm this 
observation. Van Heerden and co- 
workers postulate therefore an un- 
steady-state heat transfer mechanism 
which is dominant at the “entrance” 
section of a heater, where particles at 
bed temperature touch the hot sur- 
face. Here, the full temperature dif- 
ference between wall and bed gives 
rise to high local heat transfer co- 
efficients, a thin gas film acting as 
conductive medium. When moving 
along the wall temperature equili- 
brium between particles and wall will 
be established so that the steady-state 
conduction through the laminar film 
of the particle-fluid suspension be- 
comes more important. 


A mechanism based on _ purely 
transient heat transfer is proposed by 
Mickley and Fairbanks ‘5, namely 
unsteady-state heat diffusion into 
packets of the solid-fluid suspension 
which are in rapid exchange at the 
heater wall. 


This paper is based on experiments 
which have been undertaken to shed 
more light on the mechanism of heat 
transfer in a fluidized bed. The 
authors take the view that in most 
cases heat transfer in fluidized beds 
is made up of both steady-state and 
unsteady-state processes. This should 
show up in the heat transfer correla- 
tions through the presence of two 
additive terms. The relative impor- 
tance of the latter depends on the 
experimental conditions, e.g. the size 
and shape of the heat transfer surface, 
its position in the bed and the bed 
dynamics. It is obvious that the 
smaller the heater surface the larger 
the contribution by the unsteady- 
state process. Thus, by making the 
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—— PLEXIGLASS 


; ~~ BRASS 


COPPER 


TUNGSTEN 


Figure 1—Hot wire probe. 


heater very small, e.g. by using a hot 
wire probe as in the present case, 
the stationary term can be neglected. 
This results in a simplified correla- 
tion formula as only transient heat 
transfer has to be considered. 


Experimental 


The hot wire probe is shown in 
Figure 1. A tungsten wire 1.79 in. 
long and with a diameter of 0.0052 
in. served as a filament. Tungsten was 
chosen because its mechanical pro- 
perties are superior to those of 
platinum. The ends of the No. 10 
gauge copper support wires were 
slotted and the tungsten wire at- 
tached by pinching and silver solder- 
ing. The two copper wires were hard 
soldered to two parallel brass rods, 
0.5 x 0.25 in., 5 ft. long, separated by 
a strip of plexiglass, 0.125 in. thick. 
Two heavy stranded copper wires 
connected the brass rods to the bridge 
circuit shown in Figure 2. The hot 
wire probe formed one arm of the 
bridge. The current through the hot 
wire could be varied from about 0.5 
to 5 amperes by means of two 
variable resistors. 


The temperature-resistance _ rela- 
tionship of the probe was found by 
calibration. To that purpose the 
probe was immersed in a mineral oil 
bath with stirrer and the total resis- 
tance of filament and leads measured 
at various oil temperatures by means 
of the bridge. From these values the 
resistance of the leads, found by 
short-circuiting the tungsten wire, 
was deducted. 

A 3-in. Pyrex tube fitted with a 
stainless steel porous plate served as 
a fluidizing column. Spherical Scotch- 
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lite glass beads were used as a solid 
(Table 1). The particle diameter D, 
given in Table 1 is the arithmetical 
mean of 100 particle diameters 
measured with a microscope. 2.86 Ib. 
of solid were used for each run. 
Table 2 shows some of the properties 
of the eight fluidizing gases. 

In a preliminary experiment heat 
transfer coefficients from the wire to 
the bed were measured with the 
probe at various distances above the 
porous plate. The vertical heat trans- 
fer profile so obtained showed a flat 
maximum at about half the height of 
the fluidized bed. Therefore, the hot 
wire was held for all experiments in 
a horizontal position at a distance of 
4.7 in. from the bottom of the 
column. 


To obtain a reading at a certain 
gas flow rate the bridge was set to a 
resistance corresponding to the de- 
sired wire temperature and the cur- 
rent through the wire adjusted to 
balance the bridge. For every gas- 
solid combination the flow rates were 
chosen to cover the whole range from 
still bed to vigorously bubbling bed. 
Furthermore, current readings were 


p HOT WIRE 
PROBE 


12 VOLTS DC. 28a In 


Figure 2—Bridge circuit. 


taken at different wire temperatures 
ty, normally 112°F. and 212°F. For 
gases of high thermal conductivity, 
however, the maximum current 
carrying capacity of the bridge 
limited ty to 167°F. Bed height, bed 
temperature and the pressure drop 
through 0.5 ft. of the bed were re- 
corded, too. Every run was repeated 
once. 


Calculation of the Heat Transfer 
Coefficients 

The power dissipated by the hot 
Wire amounts to 


JPR, (1) 


which would lead to a heat transfer 
coefficient of 


JPR, 
ee ; od 
TDwl (tw — tp) 


However, not all the heat is given off 
to the bed at the surface of the hot 
wire. There is a heat loss by conduc- 
tion through the tungsten wire to the 
copper leads, which are assumed to 
be at bed temperature. 


It can be shown that the relative 
error amounts to 


hi'-h Ry, 2 
—— =~ (2) 
h} Pee, Bl 


where 


h 
io (3) 


ke De 

In Figure 3 heat transfer coefficients 
obtained with particles of a diameter 
D, = 11.5 x 10% in. are plotted versus 
G/Gry¢- The results for four gases are 
shown only. The similarity between 
the curves suggests the use of the 
maximum of every curve as a re- 
ference state to compare the thermal 
properties of the various solid-fluid 
combinations. h,,,, - values as well as 
the corresponding reduced mass ve- 
locities are found in Table 3. 


TABLE 1 
PROPERTIES OF GLASS BEADS* 


Sieve 


Manufacturer's 
Limits 


10-3 ins. 


Pb 
Ibs./(ft.3) 


ps 
Ibs./(ft.3) | 





— 150 + 200 


— 80 + 100 


109 





* 3 m “Superbrite” 


154 93.0 


152 
152 


156 


glass beanie a by Minnesota Mining and Manufacturing 


Co. Ltd., 900 Fauquier Ave., St. Paul 6, Minn. 
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Heat Transfer Mechanism 


The solids flow pattern in a 
fluidized bed is such that the particles 
can be assumed to move essentially in 
a direction perpendicular to the axis 
of a horizontal wire submerged in the 
center of the bed. In the model pro- 
posed here heat transfer from the hot 
wire to the fluidized bed is effected 
by the transient heating of particles 
touching the wire for a short time or 
passing it very closely. After having 
left the wire the particles give off 
their heat to other particles and the 
fluidizing gas. Obviously, the time 
between two successive contacts of 
one particle with the wire is large 
compared with the relaxation time of 
the particle-gas system. Consequently, 
particles arriving at the wire are at 
bed temperature. 

The following assumptions have 
been made: 

1. Heat transfer from the wire to 
the particles occurs only by con- 
duction through a gas film, solid- 
solid heat exchange being negli- 
gible. 

. Thermal contact between one 
particle and the wire is esta- 
blished during a certain time 7, 
which is the residence time of 
the particle in the vicinity of the 
wire. 7 is the average value and 
applies to all particles of a par- 
ticular fluidized bed. 

. Heat transfer by forced convec- 
tion to the fluidizing gas is 
neglected. 

. No temperature gradient is es- 
tablished within the _ particle 
during the heat-up period. 

No temperature gradient exists 
in the filament in a direction 
perpendicular to the wire axis. 


The heat flowing across the fluid 
space between a single particle and 
the heated wire in time dt is 


ki/b -A. (tw — tp) (4) 


Heat accumulated by the particle in 
time dt is 

Coan. . dt, ‘at (S) 
Equating (4) and (5) and rearranging 
ky .A. 
bCpare. ; 
Let a = k,A/bC,,,,. From an integra- 
tion over the average residence time 


r, and the subtraction of (t, — t,) 
from each side, 


dtp/a = (tw — tp) (6) 


tp — tp = (tw — th) (1 —e-@7) = (7) 
During the residence time +, the heat 
quantity received by one particle is 
q = Coart. (tp — to) = 

Chart. (tw — ty) (1 — e487) (8) 


h btu. /ft® /hr. /deg. F. 


If n is the number of particles touch- 
ing unit wire surface in unit time, the 
expression for heat transfer coeffi- 
cient can be written 


h = 


TABLE 2 


PROPERTIES OF GASES 


ks Temperature Prandtl 
at 32°F. coefficient No. 
Btu/(ft.) (hr.) of kg Com 
(deg. F.) (deg. F.)~! = 
f 


0.101 0.00155 0.718 


0.0849 0.00117 0.638 


.0719 0.00155 444 
.00155 
.0294 .00158 


.00158 


.00161 
.00261 


.00372 


©——o— Hp 149° F. 
rae 
ntti He 167° F. 


-3 
Dp = I1.5x IO ins. 


9° 
x 


ei 38% Hp 62% No I67°F 


x 


oe | 
a ~ 
a’ 
aa adi ~ ——o AIR 212° F, 


© 
nan 


Figure 3—Heat transfer coefficients vs reduced mass velocity. 


particular size group has been chosen. 
Substituting, 


n. Coase. = ho (1 — €) (10) 
TA 


: —Ak and p = — (11) 
_ nae. Cone = CS ") (9) bC part. 
=a 


bC part i 
equation (9) now reads 


Since n is proportional to the num- 
bers of particles per unit volume, it is 
also proportional to (1 — ¢«), when a 
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h = ho (1 — €) (1 — e- PK) (12) 


h, and p are experimental constants, 
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TABLE 
CALCULATED AND EXPERIMENTAL hinax - 


Dp 


3 
- VALUES 


Rimexs exp. Rimax, calc. 


Be Deviation 


tw ty ; 
10-3 ins. Gas deg. F.| deg. F. |G/Gmtex,| 1—€ Y 


| Btu/(ft.2) (hr.) ( 
| (deg. F.) | 


H, | | 118 | 0.481 | 1,260 
He | | 248 | 0.406 | 1,065 


| 81.0 He | 
19.0 No 


| 78.9 He 
| 21.1 Ne 


| 33.3 He 
66.7 Ne 


CO: 


H. 
He 


81.0 He 
19.0 Ne 


78.9 He 
21.1 Ne 


33.3 He 
66.7 Ne 


Air 
CO, 
Freon-12 | 


H, 
He 
81.0 H: 


9 
» 


19.0 N: 


78.9 He | 
21.1 Ne 


| 33.3 He 
66.7 No | 


Air | 212 
co, | 212 
Freon-12| 212 











H, | 149 
He | 167 


81.0 He 
19.0 No 


78.9 Hy 
| 21.1 Ns 


| 37.6 He 
62.4 No 


Air | 212 | 68 
CO, | 212 68 
Freon-12| 212 65 


167 | 63 
167 | 61 


167 57 





258) 


Air | | | 188 | 0.415] 513 
80) | 0.425 | 375 
| Freon-12 | ge. 145") | 0.42 275 


3 535 | 147 
15} 0.543 | 102 


230) | 0.437 | 1,115 


216 | 0.440 | 1,050 


740 


,050 
990 


960 


950 


595 


395 
305 
217 


865 
810 
840 


775 


480 
290 
230 
150 
705 
645 
650 








630 


315 
210 








(1) as calc sae according to Graf E., Thesis. Zurich (1955). 


their values depending on experimen- 
tal conditions. Using h,,,, values de- 
fined earlier, a set of constants has 
been computed for each specified 
particle size. These sets have been 
calculated from the experimental re- 
sults by the method of least squares; 
and are given in Table 4. h,,,, values 
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calculated with these constants are 
compared with experimentally deter- 
mined h,,,, Values in Table 3. 

The temperature t, at which the 
thermal conductivity k, of a gas was 
taken, was chosen to be the arith- 
metical mean between the wire tem- 
perature t, and a mean particle 
temperature t, 


TABLE 4 
SUMMARY OF EMPIRICAL CONSTANTS 
IN FORMULA (12) 


Dp | ho | p 
10-3 ins. | btu/(ft.2) ) (hr.) | (ft.) (hr) 
i 


(deg. F.) | (deg. F VY ‘Btu 


36 


| 2,660 | 

| 2,110 30 
| 

| 


1,660 25 


1,550 18 


oe 7 - mei — e- Pkt) 
(13) 

To calculate t, from equation (13) 
an approximate value of p was deter- 
mined by taking ky at the mean tem- 


perature between wire and bed. 


In Figure 4 experimental hy,.— 
values have been plotted against the 
gas thermal conductivity. The curves 
were calculated from equation (12) 
using the values for h, and p given in 
Table 4. 


Discussion 


The fact that heat transfer co- 
efficients from a hot wire to a fluid- 
ized bed can be reasonably well pre- 
dicted for eight different fluidizing 
gases with the use of only two ex- 
perimental constants is strong support 
for the mechanism proposed here. So 
far, the two constants h, and P have 
to be determined by experiment. 
Further experimental work has to be 
done to arrive at correlations for h, 
and p with variables like G, D,, D, 
and the physical properties of gas and 
solid, as well as the dimensions of the 
column. 

A more general correlation involv- 
ing the constants h, and p, should 
yield useful information about the 
dy namics of a fluidized bed. h, has 
the dimension of a heat transfer co- 
efficient. The product h, (1—e) 
actually represents the inherent heat 
transfer potential of a fluidized bed 
in contact with a hot wire. Because 
of the finite values of gas thermal 
conductivities only a fraction of this 
potential can be utilized, this fraction 
being represented by the factor 
(1—e> pk), According to equation 
(10) h, (1—e) is closely related to 
the solids mass flow past the wire: 


ho (1 — €) T 
nN Mpart. = — =n X 6° D,? ps 


Cs 
(14) 
(1— 


 & 
The dimension of —— = is that of 
8 
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Figure 4—Maximum heat transfer coefficients vs gas thermal conductivity. 


a mass velocity, ML—*t —!. There is 
not necessarily a simple relationship 
between this mass velocity and the 
solids motion of the bed, because of 
the possible interference to flow of 
particles by the hot wire itself. But 
since the solids mass velocity is re- 
lated to the linear particle velocity it 
should be possible to measure particle 


velocities with a suitable hot wire 
probe. 


Nomenclature 


A Area of heat path 

b Length of heat path 

Cyart. Heat capacity of one particle 
Ce Heat capacity of solid 

Dp Particle diameter 

Dw Wire diameter 
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Superficial mass velocity of fluid 
Superficial mass velocity at incip- 
ient fluidization 
Heat transfer coefficient 
Electric current through wire probe 
Thermal equivalent of electrical 
energy 
Thermal conductivity of fluid 
Thermal conductivity of wire 
material (tungsten ) 
Length of hot wire 

Mpart. Mass of one particle 

n Number of particles per unit time, 
touching unit wire surface 

Rep Resistance of wire at bed tempera- 
ture 

Riw Resistance of hot wire 

tp Bed temperature 
Temperature at which thermal con- 
ductivity of fluids were taken 
Particle temperature 
Wire temperature 
Time 
Void fraction 
Viscosity of fluid 
Bulk density of solid 

ps Density of solid 
Residence time 
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Counter Diffusion of Ions in Water’ 


E. R. GILLILAND?, R. F. BADDOUR®, and D. J. GOLDSTEIN* 


In most ion diffusion processes, more than one ion has a concentra- 


tion gradient. 


potential gradient by the diffusion process itself. 


Differences in ionic mobilities result in the buildup of a 


Thus, for ionic diffu- 


sion problems, the flux must be expressed as the sum of two terms, one 
involving concentration gradient and the other involving the electric 
field. Such an expression has long been known to apply to salt diffusion. 
It is shown here that the same expression is applicable to the simultaneous 
diffusion of two ions of the same charge in opposite directions, whether 
or not the third, common ion is mobile. 


The general equations for the simultaneous diffusion of many ions 


are set up. 


These equations are integrated for the case of diffusion in 


diaphragm cells. The results of this integration agree well with experi- 
mental data obtained for the counter diffusion of hydrogen ion with 
ammonium, calcium, and lanthanum ions, all in tenth normal chloride 


solutions. 


It is concluded that rates of counter diffusion of simple ions in dilute 
water solutions can be estimated quite well from the ionic equivalent 


conductances alone. 


y counter diffusion is meant the 

diffusion of two ions of the same 
charge in opposite direction. This 
problem occurs in many cases of mass 
transfer where diffusion of ionic 
species is involved. Counter diffusion 
is involved in correlating liquid side 
resistance to the rate of ion ex- 
change “). Sherwood and Wei ©) 
have recently used the methods dis- 
cussed here to investigate the rate of 
neutralization of HC] by NaOH in a 
diaphragm cell. 

These two cases will be discussed 
later, with the wider aspects of the 
problem, when a specific case has first 
been dealt with. This specific case is 
counter diffusion in diaphragm cells 
and the theoretical prediction of the 
experimental results. 


Experimental Equipment 
Experiments were made in dia- 
phragm cells. A cell is pictured in 
Figure 1. It is a device for obtaining 
two well stirred reservoirs separated 
by a completely stationary water 
layer. The two reservoirs are filled 
with different solutions and _inter- 
change between the reservoirs takes 


place by diffusion through the sta- 
tionary layer separating them. The 
concentration of the various species 
in each reservoir is determined as a 
function of time. The diffusion zone 
is held stationary by making it a 
water filled, porous, glass diaphragm. 
It is the repeated experience of many 
investigators that diffusion coeffi- 
cients determined in a diaphragm cell 
are the same as diffusion coefficients 
determined by other methods ‘), 


Procedure. (1) The lower reser- 
voir, the one with the stop-cock, was 
rinsed with the desired solution, and 
then the reservoir and diaphragm 
were filled with this solution. (2) 
The upper reservoir was rinsed with 
a different solution and some 50 cc. 
put into the reservoir. (3) The cell 
was left for 45 min. to establish con- 
centration gradients through the 
diaphragm. (4) The upper solution 
was poured off and a volume of fresh 
solution equal to the volume of the 
lower compartment was pipetted into 
the upper reservoir. The time was 
noted as the start of the run. The 
material balance at the end of the run 
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always showed that there was no 
noticeable change in the concentra- 
tions in the lower reservoir during 
the preliminary run. (5) The cell was 
put in a water bath at 25°C. and 
rotated about its longitudinal axis 
with this axis at about 7° to the 
vertical. This caused the glass sinker 
and float to roll around the dia- 
phragm and stir the reservoirs. (6) At 
the end of the run the reservoir solu- 
tions were drawn off and analyzed. 


Systems Studied. The counter 
diffusion of H+ with NH,+, H+ 
with Ca++, and H+ with Lat++ 
was studied. The anion was always 
Cl- and the initial total normality 
was 0.1. Table 1 gives a list of the 
solutions initially put into each reser- 
voir. 

These systems were selected be- 
cause they represent a good range of 
valency. Also, all the ions except 
Lat+++ could be readily determined 
chemically. H+ is the fastest moving 
ion, and these systems have a large 
difference between the ionic self dif- 
fusion coefficients of the counter dif- 
fusing ions. It was therefore thought 
that a significant movement of Cl~ 
might occur. 


The total normality of 0.1 is large 
enough to make chemical analysis 
simple and small enough so that 
infinite dilution values of diffusivities 
might be applicable. The initial nor- 
mality difference between the reser- 
voirs of 0.04 was large enough to 
prevent errors caused by differencing 
two large numbers, and small enough 
so that some idea of the effect of the 
ratio of concentrations of the counter 
diffusing ions might be obtained. The 
initial normality difference of 0.1 is 
the worst case from the point of view 
of the approximations in the theo- 
retical treatment. 


Chemical Analysis. The primary 
standards were prepared by dissolv- 
ing known weights of silver nitrate or 
sodium chloride in known volumes of 
water. The standards were titrated 
against each other. Reagent grade 
chemicals, without further purifica- 
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tion, were used. Stock solutions of 
exactly 0.IN HCl, NH,Cl, CaCl, and 
LaCl, were prepared. The various 
initial reservoir solutions were made 
by mixing appropriate volumes of 
these stock solutions. 

Hy drogen ion was determined by 
titrating with ammonium hy droxide 
or sodium hydroxide to an end point 
with bromophenol blue. Ammonium 
ion was determined by adding 5 cc. 
of neutral 37% formaldehy de to the 
neutral solution and determining the 
liberated hydrogen ion by titrating 
with sodium hydroxide to an end 
point with phenophthalein. Calcium 
ion was titrated with ethylene di- 
amine tetraacetic acid disodium salt 
containing a little magnesium ion to 
an end point with Eriochrome black 
T in the presence of an ammonium 
hydroxide-ammonium chloride buf- 
fer. Lanthanum ion was not deter- 
mined. 


Accuracy. The question of 
whether transfer through the dia- 
phragm occurs by other means than 
diffusion has been fully discussed 
(3, 4,5), It is not thought to have hap- 
pened in these experiments. The 
water bath was controlled at 25+ 
0.2°C. A change of 0.2°C. will not 
change the ionic diffusion coefficients 
by more than 1%. Also this tem- 
perature change will not cause a 
significant volume of solution to 
be pumped from the lower, closed 
reservoir into the diaphragm. Many 
of the runs of less than three days 
were made three times and the re- 
sults are reproducible. The accuracy 
is probably limited by the chemical 
analysis, and it is estimated that the 
concentrations were determined to 
within 0.0005 N. 


Calibration of Diaphragm Cells. 
If a Fick type equation 


N = —-D56c/6x (1) 


is integrated for a diaphragm cell 
with the coefficient assumed constant, 
the result is (3) 


in (Ac;/Ach) = — BDt (2) 


where f is the cell constant given by 


Afi 1 
fat + 3 
B i ( + 3) (3) 


The method of integration is similar 
to the one used later with the more 
complex differential equation for 
counter diffusion. 6B was determined 
by allowing 0.IN KCl to diffuse into 
initially pure water until the reser- 
voir concentrations were about 0.088 
N and 0.012N (48 hours). A value 
of D, the integral salt diffusion co- 
efficient of potassium chloride, of 
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Figure 1—Glass diaphragm diffusion 
cell. (Material is pyrex; all dimen- 
sions are approximate ). 


1.87 x 10-5cm?/sec (as recommended 
by Stokes (5) was used to compute 6. 


Two cells were used in this series 
of experiments. They were calibrated 
repeatedly over the course of nine 
months. The cell constants deter- 
mined from these calibrations did not 
vary over this period. Values of B 
for the two cells were 0.0844 + .008 
and 0.0842 + .008 (the error figures 
are average deviations of 6 different 
calibrations, each calibration being 
the average of three runs). A value 
of 6 of 0.0843 was used in all the cal- 
culations. 


Experimental Results 

The experimental results are pre- 
sented graphically in Figures 2 to 17. 
The legends for the graphs are 
lengthy and have been listed in Table 
| instead of on the graphs. The dif- 
ference in reservoir concentrations 
for each ion is plotted against the 
time. The results show some interest- 
ing features. The plots for the 
counter diffusing ions are nearly 
straight lines on a semi- -logarithmic 
scale. As equations (1) and (2) show, 
this means that a constant coefficient 
can be assigned to each ion satisfac- 
torily. However, the curves are not 
quite straight, being slightly concave 
upwards, showing that the apparent 
coefficient decreases as time increases, 
ie., as Ac decreases. The Lat+++ 
points were computed by assuming 
electroneutrality in each reservoir 
since Lat+++ was not determined 
chemically. 

The chloride ion, which is referred 
to as the common ion, has zero con- 
centration gradient at zero time, and 
must have zero concentration gradient 
at equilibrium at infinite time. Thus 
the concentration difference should 
pass through a maximum. This is the 


case, and the peak value of Ac is 
quite large, in some cases more than 
25% of the initial Ac for the counter 
diffusing ions. 


Theoretical Treatment of the 
Simultaneous Diffusion of 
Many Ions 


The basic theory is well esta- 
blished (3.5). It will be mentioned 
here because the underlying assump- 
tions need emphasizing. First, a clear 
distinction must be made between 
the practical coefficient, D*, and the 
ionic self diffusion coefficient at zero 
concentration, D. The practical co- 
efficient is in general a function of 
the mobilities, concentrations and 
concentration gradients of all the ions 
present. Its usefulness depends on 
how constant it is. It is demonstrably 
useful for the experiments done here, 
and the method of calculating it will 
be discussed. D is a property of the 
individual ion. It is computed from 
the ionic equivalent conductance, A, 
by the relationship 

RTA 
BY 
D is related to the ionic electrical 
mobility, u, by the equation 

pene (5) 
Fz 

It should be noted that a table of 
ionic equivalent conductances at 
infinite dilution is the only informa- 
tion needed for the calculations in 
this work. However, activity data are 
helpful, especially if a higher total 
normality is used (the experimental 
total normality was 0.1). 


D (4) 


Since all ions do not move at the 
same rate under the same concentra- 
tion gradients, most diffusion experi- 
ments result in a separation of ions, 
and an electric potential is set up by 
the diffusion process itself. Thus, in 
the general case, the equation N = 
— Déc/dx is not correct; an extra 
term involving the potential gradient 
must be added. The most useful way 
of expressing ionic flux in free dif- 
fusion is, therefore, to write. 

Ni = E ac = D; 6 Ci 6 xX (6) 
for ions of the i type. The quantities 
A, u, z, have been taken positive for 
cations and negative for anions. 
Equation (4) shows that with this 
convention, D is always positive. 
This convention means that equation 
(6) applies to cations or anions. 

Now, although a separation of ions 
occurs, it is very small. If n equiva- 


*The use of two different symbols, 
D and D, is not usual, and has been 
adopted here for clarity. 


Il 


lents of a cation are placed r centi- 
meters away from n equivalents of an 
anion in water, the potential gradient 
between them is about 1.1 x 1015 
n/r? volts/em. Thus only a very 
small separation of ions is ‘needed to 
produce a large potential gradient, 
and, in making a material balance it 
is a very good approximation to as- 
sume electro-neutrality and write 

Lei +) = Lei( —) (7) 
where © means summed for all the 
ions, cations or anions, present. Dif- 
ferentiating, 

2Gi(+) = 2Gi(—) (8) 
where G has been used as an abbre- 


viation for §c/8x. To preserve elec- 
troneutrality it is necessary to put 


ZNi(+) = ZNi(—) (9) 
Combining equation (6) with equa- 
tion (9) to eliminate E, and with 
equation (5) to eliminate u, gives 
Ni = — DiGi + 

SDiGK+) — ZDiGd-) 
=D; zici(+) — - DDizici( —) 

Vinograd and McBain ‘*) derived 
equation (10) and have successfully 
used it to predict the simultaneous 
diffusion of two salts. Sherwood and 
Wei (2) have used it to predict the 


rate of neutralization of HCl by 
NaOH and of acetic acid by NaOH. 


Diz (10) 


Integration of the Diffusion 
Equations for Counter Diffusion 
in a Diaphragm Cell 

In the simplest case of counter dif- 
fusion, which includes the experi- 
ments described here, three ions are 
involved. Let suffixes 1 and 2 refer to 
counter diffusing cations, and suffix 
3 refer to the common anion. Equa- 
tions (7) and (8) can be used to 
eliminate c, and G, from equations 
(10) to give 





N, = fGi+gGp» (11a) 
N: = hG, + kG. (11b) 
where 
cs Di D; (z3—2z1)e1 — Dy (Dez2— D3z3)c2 
(Dizi — Dszs)er + (Dez2— D323) Co 
(12a) 
_ D,; (D2—Ds)zici 
- (Dyz1 — Dzs) )ex + (Dez2— Dszs)eco 
(12b) 
De (Di —Ds)z2c2 
oS io Dae + ace Dszs)c2 
(12c) 
i D2Ds (23 — 22)c2 — De (Diz: — Dszs)er 


(Dyz1— Dszs)e1 + (Deze— Dyz3)co 
(12d) 
Now two further approximations 
will be made. The concentration 
gradients will be taken as linear, i.e., 
the substitutions 


12 


G, = Ac,/I, 
will be made, and the quantities, f, g, 
h, k, will be taken as constants during 
the integration. f, g, h, k, will be 
evaluated at average values of c, and 
c,. These values were the arithmetic 
average of the reservoir concentra- 
tions. Because the reservoirs have 
equal volume the arithmetic average 
concentration is also the concentra- 
tion everywhere at infinite time. 


G, = Ac:/1 (13) 


At the boundaries of the diaphragm, 
the concentration, c, must equal the 
reservoir concentrations, c’, c”. Also, 
in the reservoirs 


AN/V’ = 6c’/dt and AN/V” = — 6c”/ét 
(14) 

so that 

6 (c’—c”)/dt = dAc/dt = 


(A/V' + A/V”)N (15) 


Combining the boundary condition 
(15) with equations (11) and the 
approximations (13), gives 


6Aci/dt = B(fAc, + gAcez) (16a) 
bAce Ot = BihAc, + kAcz) (16b) 


The integration of equations (16) 
for the initial condition 


Ac: = — Acz = Apatt =0 (17) 
gives the result 

Ac; = A; e8 it + B, eB m! (18a) 

Acs = AgeB™t + BoeBmt  (18b) 


where m,, m, are the roots of the 
quadratic equation 

m? — (f+k)m +kf —gh =0 (19) 
and 


Ai/Ao = (m2—f+g)/(me—m;), 
B,/A. =(m,—f+g)/(mi—m_z) (20) 


= (m,—f) Ai/g, 
= (m2—f) B:/g (21) 


When these equations are applied 
to special cases, the results are reason- 
able. If c,/c, is allowed to tend to 
zero, it is found that ion 1 (which is 
present in trace quantities) moves 
with its self diffusion coefficient. The 
special case D, = 0 (which is referred 
to as the case of immobile common 
ion), and the case D, D, which 
gives the same result, will be discussed 
later. 

If D, = D,, it is found that the co- 
efficients A, ond B,, are indeterminate, 
being zero over zero. In fact, the 
solution of the equations is a little 
different. Equations (18) to (20) 
apply, but equation (21) does not 
apply and must be replaced by equa- 
tion (22 


A2/A, = h/(f—k), B2/A, = 


(k—f—h)/(f—k) (22) 


This case is important because the 
self diffusion coefficients of NH,+ 
and Cl~ are so close that equations 
(22) had to be used and not equations 
(21) in calculating the H+—NH,+ 
—Cl~ system. 


The Calculated Curves. In mak- 
ing the calculations, the following 
values of the self diffusion coefficients 
were used: “I+ = 9.311 x 10-5, 
NH, = 1.954 x 10-5, Cat + = 0.7919 
x 10-5 Latt+ = 06166: x 10-“, 
Cl— = 2.032 x 10-5, cm?/sec. These 
were computed from values of the 
ionic equivalent conductance from 
Conway’s “Electrochemical Data”. It 
must be remembered that z, is nega- 
tive. The experimental points for the 
common ion are less accurate than 
those for the counter diffusing ions, 
as the measured concentrations are 
differenced to four significant figures 
instead of three. 

Equations (18) to (22) were used 
to compute the solid lines in the 
graphs, Figures 2 to 17. The agree- 
ment between theory and experiment 
is quite satisfactory, particularly con- 
sidering the approximations involved. 
For example, for the experiments 
plotted in Figure 6, at zero time, 
c,/c, was zero at one surface of the 
diaphragm and infinite at the other 
surface, and yet a satisfactory predic- 
tion of the experimental results was 
obtained by taking c,/c, equal to one 
throughout the diaphragm for all 
time. The reason why such extreme 
variations in concentration can be 
ignored is not fully understood, but 
must be tied up w ith the fact that the 
diaphragm is very small and very 
little material accumulates in it. For 
the mathematically simpler case of 
immobile common ion, the assump- 
tion that no material accumulates in 
the diaphragm makes it possible to 
demonstrate that the variations of 
concentration can be ignored (see 
below). 


The solid lines plotted in Figures 
2, 3, 4, 5 were calculated by taking 
Cut/Cyy,+ = 0.25, 0.667, 1.5, 4, 
respectively. The agreement between 
theory and experiment improves in 
this order. Part of the reason for this 
is that the quantities f, g, h, k, change 
less rapidly with Cy+/Cyy,+ when 
Cy+/Cyy,+ is large. The same trend 
is noticeable for the Ca++—H+ 
results. 

Vinograd and McBain have suc- 
cessfully used the same approxima- 
tions as here to predict the simultan- 
eous diffusion of two salts (three ions 
moving in the same direction). How- 
ever, counter diffusion experiments 
test the usefulness of the approxima- 
tions more than salt diffusion experi- 
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TABLE 1 


SysTEMS STUDIED AND LEGENDS FOR THE GRAPHS, FIGURES 2 TO 17 


Initial Solution Normalities 


Reservoir 1 


HCl 0, NH,Cl 0.10 
HCl 0.02, NH,Cl 0.08 
HCl 0.04, NH,Cl 0.06 
HCl 0.06, NH,Cl 0.04 
HCl 0.1 

HCl 0, CaCl. 0.10 
HCl 0.02, CaCl, 0.08 
HCl 0.04, CaCl. 0.06 
HCl 0.06, CaCl. 0.04 


HCl 0, LaCl; 0.10 

HCl 0.02, LaCl; 0.08 
HCl 0.04, LaCl; 0.06 
HCI 0.06, LaCl; 0.04 


Reservoir 2 


HCl 0.04, NH,Cl 0.06 
HCl 0.06, NH,Cl 0.04 
HCl 0.08, NH,Cl 0.02 
HCl 0.10, NH,C1 0 
NH,Cl 0.1 

HCl 0.04, CaCl, 0.06 
HCI 0.06, CaCl, 0.04 
HCl 0.08, CaCl. 0.02 
HCI 0.10, CaCl 0 


HCl 0.04, LaCl; 0.06 
HCI 0.06, LaCl; 0.04 
HC] 0.08, LaCl; 0.02 
HCl 0.10, LaCl; 0 


LEGENDS FOR THE GRAPHS 


Figures 2-14: Counter Diffusing ions 


Experimental points, O = H+, @ = NH,*, Cat+ or Lat+*+ 


Calculated lines 


with immobile common ion 


Figures 15-17: Common ion 


Experimental points for Cl 
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; calculated lines 


Sop NORMALITIES 


Results Presented 


in 


Figure Number 


, and 15 curve 1 
, and 15 curve 2 
, and 15 curve 3 
, and 15 curve 4 
and 16 

, and 17 curve 1 
, and 17 curve 2 
9, and 17 curve 3 
10, and 17 curve 4 


COND Ue Wy | 


11 
12 
13 
14 


with mobile common ion and activity correction 


Figure 5 — see Table 1 


ments because the ratios of the ionic 
concentrations vary much more in 
counter diffusion experiments. 


The Effect of Finite Concentra- 
tion. The values of the ionic self 
diffusion coefficients, used to com- 
pute the solid lines of the graphs, 
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Figure 6 — see Table 1 


from 


were themselves computed 
values of the ionic equivalent con- 
ductance at zero concentration by 
using equation (4). In fact, zero con- 
centration has been assumed through- 
out this treatment. The basic flux 
equation (6) has no terms for inter- 
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Figure 8 — see Table 1 


action between the ions (other than 
the general field term) which would 
be needed at finite concentrations. 
Hartley ‘) points out that the 
chemical potential gradient is the 
driving force for a diffusion process. 
In the flux equation (6), the con- 
centration gradient term should have 
involved activities and not concen- 
trations. This means that one way to 
correct for finite concentration is to 
multiply every self diffusion coeffi- 
cient used by the correction term 
(1 — 8lny/dSlInc). Here y is the 
ionic activity coefficient, and, for the 
counter diffusing ions the variation of 
y with c in solutions of constant 
normality is needed. To get a rough 
picture of this correction, the simplest 
method was used. It was assumed, 
after Lewis and Randall, that the 
activity coefficients were functions of 
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Figure 9 — see Table | 
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the ionic strength only. This results 
in no correction for the H+ — NH,+ 
- Cl- system because the ionic 
strength is independent of C,+/ 
Cyy,t- For two sets of the H+ — 
Ca++ — Cl~ system, the correction 
was appreciable. The corrected 
curves are plotted as dotted lines in 
Figures 7 and 8. The dotted lines de- 
scribe the experiments better than the 
solid lines. In Figure 9 and 10 the cor- 
rection is negligible. Corrections for 
the H+ Lat+ + + Cl- system 
would be similar to those for the 
H+ Ca++ Cl- system, only 
larger. However, the activity coefh- 
cient of La+++ is so different from 
one that it is doubtful if the Lewis 
and Randall limiting rule can be ap- 
plied and the corrections were not 
made. 


The integration of equations (11) 
is specific for the geometry of the 
diaphragm cells. It is worth while 
considering some other approaches 
to this problem, which could be used 
for other problems. The success of 
the approximations used in integrat- 
ing suggests that the same approxima- 
tions could be used in other methods 
of solution. 


Ionic Practical Counter Diffusion 
Coefficients with Mobile Common 
Ion 

Following the method of Vinograd 
and McBain, practical coefficients can 
be defined for the two counter dif- 
fusing ions by the equations. 


Ni = «- D,Gi, Ne = - DG: (23) 


Equations (11) and (13) show that 
the practical coefficients are given by 


D: = f + (gAc2/ Aci) 
D: = k + (hAc;/Acz) 


In Table 2, under the heading mobile 
common ion, some values of D, and 
‘D, are given for two of the systems 
studied. These values were computed 
by taking the value of Ac,/Ac, at 

the very beginning of the experiment, 
at zero time; ie., Ac,/Ac, = —l. 
Also, as before, f, g, h, k were com- 
puted by using average values of c, 
and c,. If these coefficients were to 
be used in the cell equation, (2), to 
find the curve of Ac against time, the 
curves for the counter diffusing ions 
would coincide with the solid lines in 
Figures 2 to 14 at zero time and be 
straight lines having the same slope 
that the solid lines have at zero time. 
Because the experimental lines are 
remarkably straight, the constant 
practical coefficients are very useful, 
and certainly their use is simpler than 
using the integrated equations, (18) 


(24a) 
(24b) 


> 


Oce NORMALITIES 


TIME (Hours) 
Figure 16 — see Table 1 


80 


to (22). For short experiments, this 
approach to predicting the curves is 
probably the best. 


Counter Diffusion Coefficient with 
Immobile Common Ion 

If an experiment is imagined in 
which the Cl~ ion is replaced by an 
immobile ion, then only two ions are 
moving and the equations are mathe- 
matically much simplified and are in- 
teresting to study. This system is of 
some practical importance for cor- 
relating rate coefficients in ion ex- 
change ‘!), In this case there are two 
exchanging, counter diffusing ions, 
but the common ion of the solution 
is excluded from the resin, and it is 
probably best to assume it to be 
immobile. 

The two counter diffusing ions 
have equal and opposite initial con- 
centration gradients and must there- 
fore have the same diffusion co- 
efficient if electroneutrality is to be 
preserved. If, in equations (12), D 
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TIME (HOURS) 
Figure 17 — see Table 1 


is put equal to zero, (or if D, is put 
equal to D,, for in this case no elec- 
tric potential can be set up by the 
diffusion process and the common ion 
should not move), the single counter 
diffusion coefficient is found to be 
given by 

Z1C1/(Z1C1 +2Z2C2) 

Dz 
Z2C2/(Z1C1 +Z2C2) 
D; 


5, Dry _ 
(25) 


The same result can be obtained more 
directly by writing equation (10) for 
ion 1 and eliminating c, and G, by 
using the electroneutrality conditions 
(7) and (8). 

The dashed lines in Figure 2 to 14 
have been computed from equation 
(25) together with the cell equation 
(2). The lines are slightly curved 
because the dependence of 3D,, 
c,/c, has been taken into account as 
discussed below. It will be seen that 
the mobility of the common ion al- 
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TABLE 2 


CALCULATED PRACTICAL DIFFUSION COEFFICIENTS 


The Initial Reservoir Solutions were 0.1N in Cl- 


Initial Reservoir Normalities 


of H+ 
Reservoir 1 Reservoir 2 
0 0.04 
0.02 0.06 
0.04 0.08 
0.06 0.10 
0 0.10 


ways increases the apparent coeffi- 
cient of the faster of the counter dif- 
fusing ions, and decreases the ap- 
parent coefficient of the slower. This 
must be because the faster ion takes 
some of the common ion with it. The 
fact that equation (25) is obtained 
from the integrated equations which 
fit the experimental results, coupled 
with the fact that the experimental 
results always lie to either side of the 
curve calculated from equation (25) 
suggests that this equation would be 
correct if a truly immobile common 
ion had been used. 


It was mentioned that this case is 
simple, mathematically, because only 
two ions are involved, and that some 
allowance for the variation of D with 
c can be made for this case. The 
method used is Stokes’ (5). The as- 
sumption of a linear concentration 
gradient is still made, and the addi- 
tional (most reasonable) assumption 
that the diaphragm is too small for 
any material to accumulate in it, is 
made. With these assumptions, Fick’s 
equation can be integrated for the 
geometry of a cell, taking into ac- 
count that c,/c, and therefore D,, 
vary both with position in the 
diaphragm, and with time. 

Stokes’ treatment is as follows. If 
no material accumulates in_ the 
diaphragm, the flux must be constant 
through the diaphragm. That is, 
‘Ddc/dx is independent of x, and a 
coefficient ‘D’ may be defined by the 
equation 


D’ Ac/l = Déc/6x (26) 
‘(—D' is a function of time, and there- 
for of Ac, but is independent of x. 


Now suppose that Fick’s equation, 
(1), is to be applied to a diaphragm 
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The Cations Being H+ and Either NH,*+, Ca++ or Lat+t++ 


Calculated Practical Diffusion Coefficients, cm?/sec x 105 


Mobile Common Ion (eqn. 24) Immobile Common Ion 


NH,t+ — H+ Cat+ — H+ 
H+ NH, H+ OY [ee sticteea tn = ieee 

=0 t=oa =0 t=a t=0 t= © 
6.80 4.08 | 6.24 2.89 5.68 | 5.30 4.73 4.24 4.91 4.46 
5.36 3.20 4.61 2.00 3.83 
4.42 2.64 3.54 1.44 2.90 2.85 ee 1.65 1.70 1.63 
3.76 2.20 2.84 1.05 2.34 2.32 | re | 1.14 1.06 1.03 


4.84 2.90 4.00 1.69 3.86 gee 2.76 2.03 2.81 2.06 


(ii) It is assumed that no material 
accumulates in the diaphragm. In this 
case, material conservation says that 
the mass of ion i in the two reservoirs 
is constant. Since, in all the work 
here, the two reservoirs have equal 
volume, the concentration of ion i in 
the two reservoirs may be put equal 
to a constant. 


cell, but D is a differential coefficient, 
a function of c. The substitution 
5c/8x = Ac/l can be made if D’ is 
used instead of D. When this is done 
and the boundary condition (15) is 
used, the cell differential equation is 
obtained 


6Ac/ét = — BD'Ac (27) 


The integration of this equation for 


nly ; ‘ + 2 ere, aa ? 
constant ‘D’ gives an equation like Let 2b1 = cr +e (32a) 
(2). However, D’ in fact varies with 2be = co’ + co” (32b) 
time, and so an integral coefficient : 
‘H 8 aan b,, b, are the average concentrations 
D, will be defined by the equation of the ions. Because the two reservoirs 


have equal volume, b,, b, are also the 
concentrations of the ions which 
prevail everywhere after an infinitely 
long experiment. 


In (Ac;/Acp) = — BDt (28) 


and ‘(D will be related to D’. To arrive 


at ‘D, ‘D’ is first found. From equation 
(26) 


Since equation (31) holds in the 
reservoirs as well as in the diaphragm, 


a 


‘D’ = (1/Ac) | Dde (29a) 


” 
Cc 


by + bo =a (33) 


Using equation (31), D can be ex- 


Then from equations (28) and (27) pressed as a function of c, only. 


? D = D, De [(21—z2)e1 + zeal 
‘D = in (Ac; /Acy) J caac D' Ac) = (z:1D:—zeD;)c1 + zoD oa 
Cb Then from equation (29a) 
(29b) 
: : us DD». (2, —Z2) 
In the counter diffusion case, the D' = 


differential coefficient ‘D, from equa- Z1Di—22Dz 
tion (25), can be rewritten aZ1Z2D),Dz2 (D;— D2) 


1 
A, (z1D;—z2D2)? . 


(21D, —z2D2)ce)’ + z2Dea 
(z:D; —Z2D»2)c,” + z2Doa 


D, Dz (z1¢1 +222) 


(30) 
Z1C,D, +zc2D». 


It is best at this point to introduce ] 
two new svmbols. where A has been written for c’ -c”. 

Since c,’ = b, + A,/2 and c,” = 
b, — A,/2, the log term in D’ can be 
written 


(i) Only two ions (which have the 
same charge) are moving and electro- 
neutrality insists that c, + c, be a 


constant everywhere. [2 z2Dsa/(zi\D:—z2D2)] + 2b; + A: 


1 eee eee a 
Leta =ate (31) ” [2 2sDsa/(ziD,—2:D,)| + 2b; — di 
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Define 
Q = 2z2Dea/(z1D,—z2D2) + 2b, (34a) 
P = az,\z2D,De2 (D; — De2)/(z;D1 —z2D2)? 
(34b) 
S = D,D» (2, —2z2)/(z1D1 —z2D2) (34c) 
Then 
=§ + (P/A) in(Q + A)/(O — A) 
(35) 
Incidentally, when z, = z,, S = 0, 


and in fact z disappears from the 
equations and the algebra becomes 
very much simpler. 

In general, D must be found from 
‘D' by graphical integration of equa- 
tion (29b). However, the two ex- 
treme limits of DM when a = O and 
when t = o, can be evaluated by 
taking limits. 

1. t —> 0, Ar = A, 
D=S+(P Ap) In (Q + Ap) (QO — Ab) 


(36) 
2. Whent > o, A, > 0. 


D =S+2P/Q (37) 
If, in equation (37), the definitions 
of S, P, Q are substituted from equa- 
tions (34), and a is eliminated by 
equation (33), it is found that 
Zb1 ‘(zibi + Zob2) 

= + 

De 


Zobe/(z1b; + zeb») 
D, 


1/D = 
(38) 


This is like equation (25), and 
shows that the integral coefficient for 
an infinitely long experiment is the 
same as the differential coefficient 
computed at the values assumed by 
the concentrations everywhere at 
infinite time. 


3. It can be shown that the integral 
coefficient always lies between the 
values given by equations (36) and 
(37). 


Table 2 shows some values of 
counter diffusion coefficients with 
immobile common ion for the three 
cation pairs used in the experiments. 
In calculating the dashed lines on the 
graphs, D was assumed to be a linear 
function of Ac between its values at 
Ac, (¢ = o) and Ac = 0, 
(t = o). A few graphical integra- 


NOTE: 
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tions of equation (29b) show that 
this is true to within about 15%. 
The notable feature of these results 
is that ‘D varies little with time. When 
this work was being done, it was the 
discovery that ‘D, in this case, varied 
so little with time, that encouraged 
the use of constant values of con- 
centration when _ integrating the 
equations for a mobile common ion. 


Conclusions 


The experimental results show that 
the electric potential set up by the 
counter diffusion process can cause 
the common ion to build up a con- 
siderable concentration gradient when 
it had no concentration gradient 
initially. The plots of concentration 
difference against time on a semi- 
logarithmic scale are nearly straight 
lines for the counter diffusing ions. 
This means that a single practical 
coefficient can be assigned to each of 
these ions to get a satisfactory pre- 
diction of the experiments. 

The agreement between theory 
and experiment is most satisfactory. 
It is not thought that this agreement 
is specific to the sysems studied. Any 
similarly simple system would work. 
By a simple system is meant simple 
ions in dilute solution. No allowance 
has been made for complex ion 
formation, or partial ionization, al- 
though the success of Sherwood and 
Wei in predicting the counter dif- 
fusion of acetic acid and sodium 
hydroxide suggests that partial ioni- 
zation can be allowed for. The 
activity corrections made were of the 
simplest sort and would be _ inade- 
quate for more concentrated solu- 
tions than those used. 

It has proved possible to predict 
counter diffusion results using only 
a table of ionic equivalent conduc- 
tances at infinite dilution. The 
method used is already known to 
work for salt diffusion, but counter 
diffusion is a greater test of the 
approximations in the theory. 


Nomenclature 

a total cation concentration of each 
reservoir. In the « case of immobile 
common ion, a = ¢;)’ +c.’ = ¢)” 
+ ce” =e, + C2. 

A effective cross section of the dia- 


phragm, cm? 
AAs coefficients used in the mobile 
common ion equation (18) 


2b total ion concentration in the cell. 
Because the diaphragm has negli- 


gible volume, 2b = ¢,’ 

2be = C2’ + Ce. 

coefficients used in the mobile 

common ion equation (18) 

c concentration, equiv./cc. If the 
numerical value is followed by N, 
the units are normalities, equiv./ 


ers 
Bi,Bs 


litre. 

ee concentration in the two reservoirs 

D ionic self-diffusion coefficient, at 
zero concentration, cm?/sec. 

D practical diffusion coefficient, cm? 
sec. 

D integral practical diffusion coef- 
ficient 

D' coefficient averaged through the 
diaphragm 

E potential gradient, volts/cm. 

f,g,h,k quantities defined by equations 
(12) 

F The Faraday, 96,488 coulomb/ 
equiv. 

G abbreviation for 6 c/6 x 

l effective diaphragm thickness, cm. 


m;,M2 exponential indices used in equa- 
tions (18) 
N flux, equiv./(cm?) (sec. ) 


P,Q,S quantities defined by equations 
(34) and used for the case of im- 
mobile common ion 

R The gas constant, 8.314 joules/ 
(gm —ion) (°K) 

¥ temperature, °K. 

t time, secs. 

u ionic electrical mobility, em?/(sec) 
(volt). Positive for cations, nega- 
tive for anions. 

V'.V” reservoir volumes, cm*. 

x distance through the diaphragm 
along the diffusion path, cm. 

z ion valency, equiv./gm—ion. Pos- 
itive for cations, negative for 
anions. 


6 ll t Se + 2 
‘ell constant = >m-? 
| cell constan iy \" cm 


A,Ac =c’-—-c 

nN ionic equivalent conductance, cm? 
(equiv) (ohm). Positive for cations, 
negative for anions. 


Suffixes 

b beginning, at time zero. 

f final, at time t 

i refers to ions of the i type. 
12 counter diffusing cations. 
3 common anion. 
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The Effect of Sonic Energy on 


Mass Transfer in 


Solid-Gas Contacting Operations 


J. W. HODGINS?, T. W. HOFFMAN*, and D. 


The addition of sonic vibrational energy to the gas in solid-gas 
contacting equipment causes a significant increase in the rate of transfer 
of material between these two phases. The amount by which the mass 
transfer rate is raised is a function of the gas flowrate, with greatest 
effects for low flowrates through fixed beds. 


Mass transfer results are presented for the evaporation of naphtha- 
lene from the surface of coated glass beads into a stream of air passing 
through them with the conditions varying from fixed bed to a fully 
fluidized bed. Experiments with a single naphthalene sphere are de- 
scribed, and mention is made of a hot-wire microphone developed for 
the purpose of measuring the sound. 

Sound intensities up to 160 decibels were used at frequencies up to 
12,000 cycles per second. Largest effects were observed at resonance 
with low frequencies (100-700 cycles per second), and the minimum 
effective intensity was 120 decibels. Maximum improvement of the mass 
transfer rate was about 120% for beds of spheres, and about 220% 
for the single sphere. This difference in enhancement effect is attri- 
buted to the acoustic properties of the cavity, and some elementary 
aspects of a cavity design are discussed. 


Introduction 
N furnaces where powdered fuel is 


probably never be reached by altera- 
tion of the flowrate, because slugging 


being burned, dramatic increases 
in combustion rate and efficiency 
have been reported ‘!) under the con- 
dition wherein sonic resonance is 
established in the combustion cham- 
ber. The improvement in combustion 
is attributed to the “scrubbing” off 
the gas film at particle surface and 
chamber wall, allowing ingress of 
fresh air to the fuel and enhanced 
heat transfer. The burning of pow- 
dered fuel is analogous to the situa- 
tion in equipment for contacting gas 
with a fluidized solid. In this case, 
however, there are distinct limits im- 
posed by the system to the amount 
of turbulence which can be intro- 
duced by mechanical means, and the 
condition of sonic resonance will 


or extreme elutriation will intervene. 
It seemed reasonable to assume that 
sonic energy might decrease surface 
film resistance to mass transfer, how- 
ever, if the sound were added to the 
gas from an external source, rather 
than producing it by the passage of 
gas. For an input power of 1 watt 
per sq. cm. cross sectional area, at a 
frequency of 256 cycles per second 
(middle C), air particles will vibrate 
with an amplitude of 5 mm.—a dis- 
tance larger than the diameter of 
most fluidized particles. Since signi- 
ficant gas movement could be ob- 
tained in this way, it was decided to 
extend the examination of sonic 
effects to fixed, as well as fluidized 
beds. 
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C. PEP 


Gas-solid contacting is used for 
many purposes—for example, catalytic 
reactions such as cracking of gas oil, 
roasting, drying, and adsorption of 
gases and vapours. In many of these 
operations the slow step is due to the 
resistance offered to the transfer of 
material between gas and solid, by 
the film at the interface. The thick- 
ness of this film is a function (among 
other things) of the degree of tur- 
bulence at the interface. Accordingly 
one test of the efficacy of sonic 
energy input to gas- -solid contacting 
equipment, is to measure the rate of 
transfer of material from the surface 
of the solid into the gas stream in the 
absence and in the presence of sound. 
In order that the characteristics of 
the bed of solids remained substan- 
tially unchanged during the experi- 
ments, it was decided to coat solid 
particles (glass beads) with a ma- 
terial (naphthalene) which could be 
evaporated into an air stream. This 
ensured that the size distribution, 
surface area, density, voidage and 
sound absorption characteristics re- 
mained constant over the experimen- 
tal period. Accordingly, this paper 
will deal with the improvement of 
the rate of transfer of naphthalene 
from the surface of glass spheres into 
an air stream, by means of the appli- 
cation of sonic vibrational energy to 
the air. 


The effect of sound energy upon 
the fluidizing property of normally 
non-fluent beds of solids has already 
been described by Morse ‘?), who 
found that low frequency (50-500 
cycles per second) high intensity 
(greater than 110 decibels) sonic 
energy will cause certain fine par- 
ticles to fluidize under conditions 
where they are normally non-fluent. 
No effect was observed on coarse, 
free-flowing granular materials. The 
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echnique of measuring mass transfer 
rates by naphthalene evaporation has 
also previously been described by 

u Chin Chu and co-workers (8), and 
bs Plewes et al. 4). Of the various 
suethods available for the presentation 
of mass transfer data, that used by 
ju Chin Chu et al. (3) was chosen, 
wherein the term 


Ke - PBM (Ng)? 3 
Gn 


as proposed by Chilton and Colburn 
‘5) was plotted against 


> = 


Dp P Vsp 
Me a — €) 
as employed by Carman ‘®). 


vs = 
Re 


Experimental Procedure 


(1) Mass transfer from coated 
glass beads 

Figure 1 illustrates the apparatus 
used in this investigation. Since the 
compressed air from the laboratory 
lines was always at a lower tempera- 
ture than the room, it was brought 
up to room temperature before its 
passage through the bed, by passing 
it through a heat exchanger (B) with 
room air blown outside the finned 
tubes. Since the laboratory tempera- 
ture was automatically regulated, air 
temperatures were thus held quite 
steady during the short runs. Quick- 
opening valves (D) gave instantaneous 


air flow and pressure relief, with 
close control of flow provided by a 
needle valve (A). The packed section 


of 4-in. glass pipe (E) below the 
column was used to minimize chan- 
nelling of the air flow. The bed was 
supported on a stainless steel screen 
(F), and the section (G) was re- 
movable for weighing. Air pressures 
above and below the bed of beads 
were recorded, and temperatures 
were measured by means of calibrated 
copper constantan thermo-couples 
(T,, T,). Air left the column through 
a manifold which took the form of a 
collar near the top, with the air enter- 
ing the collar by means of four 
equally spaced j-in. holes drilled 
the column al The sound was 
passed down the column by means 
of a 30-watt University transducer 
(H), which operated from a_ flat- 
response amplifier receiving a signal 
from a quartz oscillator. 

Sound intensity at the bed level was 
measured with a_ hot-wire micro- 
phone developed and calibrated for 
this specific purpose. This instrument 
measured the sound as a function of 
the increased heat transfer from a 
hot wire due to the agitation of the 
air. Its construction and calibration 
will be described in detail in a sub- 


To oscillator 
and amplifier 


Figure 1—Apparatus for measuring mass transfer in a bed of spheres. 


To oscillator 


and amplifier -—— R22 a 


To thermostat 
a 


celeste 


Figure 2—Apparatus for measuring mass transfer at the surface of a single sphere. 


sequent paper. Its distinguishing 
feature was its tiny size, enabling 
point measurements to be made with- 
out altering the physical characteris- 
tics of the reactor cavity. By measur- 
ing sound intensity over a range of 
frequencies, the resonant frequency 
could be established to an accuracy 
of about 0.3% 

The solid bed consisted of coated 
glass spheres 3-6 mm. in diameter, 
depending upon the value of N’,, de- 
sired. In order to make the naph- 
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thalene adhere to their surfaces, they 
were first etched by fluidizing them 
in a copper column with an air stream 
containing HF vapor. Molten naph- 
thalene was then sprayed on the beads 
as they tossed in a baffled cylinder. 
The spray was achieved by means of 
a glass atomizer, and naphthalene, 
atomizer and incoming air were all 
kept above the melting point of 
naphthalene. Spraying was continued 
until the coating represented about 
10% of the total weight of the coated 
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Figure 3—Effect of sound (600 ¢.p.s.) on mass transfer from a bed of spheres. 


beads (determined by Soxhlet extrac- 
tion with ether). 

For an experimental run, 70-100 
grams of coated beads were weighed 
in the pera 4-in. section of the 
column (G). The section was imme- 
diately aaee in the apparatus, and 
the sound (if applicable) and air 
turned on simultaneously. The inlet 
air temperature, flowrate, and_pres- 
sure were recorded. The outlet air 
temperature was recorded on a strip 
chart potentiometer, and the average 
temperature for the run was subse- 
quently determined by graphical in- 
tegration of the record. From the 
weight of naphthalene evaporated 
from the beads, the mass transfer co- 
efficient could be calculated. These 
results were obtained for a range of 
flowrates, in the presence, and in the 
absence of sound. Measurements were 
made to determine the effect of fre- 
quency and intensity of the added 
sonic energy on the rate of material 
transfer from the solid to the gas 
phase. 

Early in the experimental work, it 
was found that the first run after the 
beads were coated frequently yielded 
erratic results, due presumably to the 
elutriation of naphthalene dust. Ac- 
cordingly, before any quantitative 
experiment in which freshly-coated 
beads were used, the beads were 
fluidized for a couple of minutes at 
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high air flowrate and high sound in- 
tensity, to shake off any loosely-ad- 
hering particles. This precaution 
made for consistency in the results 
of subsequent duplicate runs. 


(2) Mass transfer from a single 

naphthalene sphere: 

A simpler system in which to in- 
vestigate the mass transfer in a 
sound beam, was to use a single 
sphere of naphthalene suspended in 
an air stream. Because of the small 
surface area, the experiments had to 
be carried on over a longer period of 
time, and this necessitated close con- 
trol of temperature during the runs. 
The apparatus is shown in Figure 2. 

Naphthalene spheres of the desired 
diameter were prepared by ev aporat- 
ing pressed naphthalene moth balls in 
a stream of hot air while contained in 
a rotating cylinder; this allowed 
evaporation to proceed uniformly. 
Final accurate sizing was accom- 
plished by rolling the sphere through 
a series of steel dies. The sphere (A, 
Figure 2) was suspended on a very 
fine steel wire at the axis of a vertical 
Lucite tube, 14 in. in diameter, and 
32 in. long (B). Surmounting the 
tube was a University 30-watt trans- 
ducer (C) which received a signal 
from an_ oscillator-amplifier com- 
bination. Metered air was passed 
through the jacket of a heat ex- 
changer (E) into a 6-litre cylindrical 


damping volume (F), covered with 
heavy insulation. A water jacket (G) 
surrounding the bottom section of the 
Lucite tube provided final control of 
the air temperature immediately be- 
fore it reached the naphthalene 
sphere. The water in the heat ex- 
changers was in continuous circula- 
tion from the thermostat, and_ final 
air temperature control was + 0.1°C. 
Temperature was measured at the 
sphere by a calibrated copper-con- 
stantan thermo-couple and con- 
tinuously recorded by _ strip-chart 
potentiometer. In order to have visual 
means of following the progress of 
the evaporation, the image of the 
sphere was projected on a target at 
intervals during the run, and thus 
qualitative observations could be 
made without interrupting the ex- 
periments. 


Sound levels were measured by 
the previously-mentioned hot-wire 
microphone, which was also used to 
determine the resonant frequency. 
The power input to the horn was 
read on an electronic wattmeter, and 
kept at a constant value manually. 


To perform an experiment, the 

following steps were followed: 

(i) The sphere and its hook were 
weighed on an analytical bal- 
ance. 

(ii) The sphere was attached to the 
support rod on the adapter at 
the top of the tube, and the 
sphere, adapter and horn fixed 
in place. 

(iii) The air and the signal to the 
horn were started simul- 
taneously. 

(iv) Time, atmospheric pressure, and 
tube temperature were re- 
corded. 

(v) At the end of the run, the loss 
in weight of the sphere was 
determined. 


Experimental Results 
(i) Experiments with beds of 
coated beads: 

Typical data for these experiments 
are shown in Table 1 (Page 21): 
(Runs marked with an asterisk were 
performed using sound at 600 c.p.s. 
and 145 decibels. ) 


Results for many experimental runs 
are summarized in Figure 3 and 
Figure 4, using the experimental con- 
ditions shown. For beds of beads 
is seen that maximum enhancement of 
mass transfer by sound occurs at low 
values of N're when the bed is sta- 
tionary, and is more pronounced a 
low frequency. 
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TABLE 1 


Wt. Surface Bed 


Jurati psa Ti ; zp 
Loss area Duration Voidage| Ht. Temp. | x Re 
(cm.) 


grams)|  (ft.2) (hours) °C.) 


0. ‘ .692 0.050 48 22. 250 
"@: . 584 0.033 .49 23. 260 
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Table 1—Effect of sonic vibrational energy on mass transfer from a packed bed 
of spheres. 


TABLE 2 


Sphere Air 
dia. flowrate 
(inches c.f.m. 


C 
€ 


Enhancement 


Ke (quiet run) 


(lb. moles hr.~! ft.~? atm.~') Ke (using sound) 


0.50 . 302 .97 220 
0.75 .318 .99 211 
1.0 . 366 .03 182 
1.50 .414 .98 136 
2.6 . 504 .O1 100 


| 


ot pe fh he 
dO hd DO be be 
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Table 2—Effect of sonic vibrational energy on mass transfer from a single sphere. 


At the 590 c.p.s. resonance point, frequency closed-pipe resonance also 
there are two lines which describe occurred, with the bed of solids as 
the mass transfer during sound input, the closure. It is also apparent that 
with a change of slope at N’,, = 860, the use of a transducer having a one- 
the point of incipient fluidization, In inch throat, attached directly to a 
all cases, no enhancement was ob- three-inch tube is inefficient, and in 
served at N’p, > 2500. addition makes prediction of the 
resonance frequency difficult. Ac- 


For quiet runs, the line relating 
cordingly, a 1}-in.-diameter column 


mass transfer factor to N’,, compares 
very well with the relationship pro- 
cured by Ju Chin Chu et al. (*), The 
experimental scatter of the above re- 
sults is considerably less than the 
previous work, and an altered equa- 
tion is suggested for the relation be- 
tween mass transfer and flowrate: 
jp = 1.051 (N’pre)~°35 for 4000 > N’re > 200 
For a fixed arbitrary value of 
N’p- = 460, the relationship between 
mass transfer factor and sound in- 
tensity is shown in Figure 5, (fre- 
quency = 590 c.p.s.). It has been a 
common observation that no effect is 


observable at intensity lower than 
120 decibels. 


Prompted by the observation that 
much larger mass transfer enhance- 
ment was obtained with _ single 
spheres at low resonant frequencies, 


an attempt was made to establish a @ With sound (145db) f= 2000cps 


resonance point in the tube at lower 0 Without sound 
frequency, but difficulty was en- 
countered in producing resonance at 
a frequency near the calculated value. 
It was soon realized that due to the 
dimensions of the first experimental 
apparatus two simultaneous standing 
waves were present in the tube. Part 
of the energy was sustaining open- 


was used, with the top of the bed of 
solids located 3 in. below the trans- 
ducer opening. This alteration in 
geometry made for more efficient 
coupling and made it possible to set 
up open-tube resonance while elimin- 
ating closed-pipe resonance. 

Mass transfer measurements at 148 
decibels and 340 c.p.s. are illustrated 
in Figure 6. 

It is interesting that, although the 
enhancement is pronounced at values 
of N’,. below 900, for higher values 
the relation between mass transfer 
factor and flowrate is the same as for 
the three-inch column. 

Figure 7 illustrates the relationship 
between pressure drop across the bed, 
and flowrate in the three-inch column 
for a shallow bed of 3 mm. naph- 
thalene-coated glass spheres. 


(ii) Experiments with a single 
sphere of naphthalene. 

The rate of removal of naphtha- 
lene proved to be strikingly de- 
pendent upon frequency. Although 
for the tube used, resonant ee 
cies occurred at 275, 480, 640, 1320, 
and 1610 cycles per second, signifi- 
cant enhancement of mass transfer 
rate occurred only at the lowest of 
these frequencies. The mass transfer 
coefficient in a sonic beam (fre- 
quency = 275 c.p.s.; intensity = 160 
db.) varies but little over a fourfold 
range of flowrate; the mass transfer 





tube' resonance, while at the same Figure 4—Effect of sound (2000 c.p.s.) on mass transfer from a bed of spheres. 
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Figure 5—Effect of sound intensity on mass transfer from a bed of spheres. 
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Figure 6—Effect of sound (340 c.p.s.) on mass transfer from a bed of spheres 


in a tube from which closed-pipe resonance has been eliminated. 


for corresponding quiet runs, of 
course, increases. The values are 


22 


given in Table 2. 


A typical trace of the projected 


profile of a single naphthalene sphere 
is shown in Figure 8. The effect of 
the sonic vibrational energy was in- 
variably that of making the evapora- 
tion from the surface of the sphere 
uniform, rather than having it con- 
centrated on the leading and trailing 
surfaces. Maximum enhancement oc- 
curred when the sphere was located 
at an antinode in the standing wave, 
and no improvement was obtained at 
the position of minimum velocity 
change—at the node. 

The results of an experiment to 
determine the effect of sphere di- 
ameter upon the mass transfer co- 
efficient in the presence of sound, are 
shown in Table 3. 


Discussion 

Effects other than increased turbu- 
lence at the gas-solid interface could 
cause an improvement in the mass 
transfer rate. For example, the values 
of the diffusivity and of the partial 
pressure of naphthalene were calcu- 
lated on the assumption that the 
temperature of the air effluent from 
the bed represented the bead surface 
temperature. In all cases where the 
sound beam was used, the effluent air 
did show a temperature rise (0.1 to 
0.4°F.), indicating some energy dis- 
sipation at the bed. At low flowrate 
one might suspect that the heat trans- 
fer coefficient at the surface of the 
beads was sufficiently low to allow the 
bead surface temperature to rise 
appreciably above that of the effluent 
air. That this was not the case was 
demonstrated by measuring surface 
temperatures during experiments, 
using fine thermocouples attached 
directly to beads in the bed. 


In every case, the temperature of 
the effluent air represented the bead 
surface temperature within 0.1°F. 
The same was true of the single naph- 
thalene spheres. 


The total pressure in the column 
increases with flowrate. Thus the de- 
crease in enhancement at high flow- 
rate could be caused by the less 
effective operation of the transducer 
diaphragm. By repeating some of the 
low flowrate experiments at elevated 
pressures, this effect was shown to be 
negligible over the pressure range 
involved. 

The increased pressure drop across 
the fixed bed in the presence of 
sound (Figure 7) is interpreted as a 
manifestation of increased turbulence 
in the air passing through the bed, 
with the maximum pressure change 
occurring at the low flowrate. 
Fluidization of course occurs at a 
lower flowrate in the presence of 
sound, and once the bed has become 
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fully fluidized, the pressure drop is 
not affected by sound input. This be- 
haviour is consistent with the mass 
transfer measurements. The reason 
that there appears to be some slight 
mass transfer enhancement in the 
fluidized bed is the method of plot- 
ting the data; the slightly increased 
voidage due to the sound has resulted 
in a higher value for N’,, for a given 
value of the mass transfer factor. 
However maximum effect was ob- 
served at a flowrate sufficiently low 
that the voidage was unaffected by 
the sound. When fluidization de- 
velops, increased turbulence mani- 
fests itself as kinetic energy used up 
in expanding the bed. The velocity 
gradients at the particle surface 
change but little, thus affecting mass 
transfer at the boundary layer only 
slightly. 

Maisel and Sherwood ‘7: 8) discuss 
the effect of scale and intensity of 
turbulence on mass transfer in tur- 
bulent air streams. In their experi- 
ments turbulence was introduced by 
placing drilled plates in the air duct 
at various distances upstream of solids 
from which the rate of evaporation 
of water was measured. Scale of 
turbulence was measured by a stan- 
dard technique involving two 
separated wire probes, and is a num- 
ber having the dimensions of length, 
which is indicative of the size of the 
eddies in the air stream. Intensity, a 
measure of the velocity fluctuations 
at a given point, was measured with 
an extremely delicate hot wire 
anemometer providing a signal to an 
oscilloscope. In all of their experi- 
ments (in which the superficial flow 
was turbulent), the effect of turbu- 
lence intensity was marked, but that 
of scale was of little consequence. 

In air in which sonic resonance is 
being maintained, the amplitude of 
the movement of the air particles is 
inversely proportional to the fre- 
quency for a given power input, but 
the displacement velocity is inde- 
pendent of frequency. Since the mass 
transfer enhancement is so very de- 
pendent upon frequency, it is ap- 
parent that the amplitude of move- 
ment of the air by the sound has a 
much greater effect than the velocity 
of the air. Interpretation of these 
effects in the light of normal turbu- 
lence phenomena such as examined by 
Maisel and Sherwood is questionable, 
for the velocity imparted to the air 
particles by the sound is 10-20 times 
greater than the superficial air ve- 
locity through the tube. An increase 
in the superficial air flowrate in the 
tube did not cause any change in the 
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Figure 7—Effect of sound on pressure drop across a bed of 3 mm. diameter glass 


spheres. 
TABLE 3 
Air velocity: 1.4 ft. per second 
Sound frequency: 268 c.p.s. 
“intensity: 163 decibels 
Dia. of sphere Ke Ne (’ ) 
(inches) (lb. moles hr.~! ft.~? atm.~') Me 
0.297 1.26 203 
0.359 1.31 246 
0.390 1.24 268 
0.453 1.11 314 
0.484 1.08 335 
0.574 0.895 410 


Table 3—Effect of sphere diameter on mass transfer at constant sonic input. 
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mass transfer rate at constant sonic 
power. 

The mass transfer coefficient de- 
creases with increasing sphere dia- 
meter for constant sound input and 
superficial air flowrate (Figure 3). 
This is to be expected for turbulent 
flow past the spherical surface—the 
mass transfer coefficient should be a 
function of (D,)*. 


Mass transfer enhancement is not 
so sharply dependent upon frequency 
when measurements are made on a 
packed bed rather than on a single 
sphere. The reason for this difference 
in behaviour is not apparent. Further 
work on these effects should be ac- 
companied by measurements of local 
turbulence effects. 


Conclusion 

Addition of sonic vibrational 
energy at low frequency and at in- 
tensity greater than 120 decibels sig- 
nificantly increases the rate at which 
material is transferred across solid- 
gas interfaces in a bed of solids. The 
enhancement is greatest for low flows 
of gas through fixed beds, and the 
improvement is negligible when full 
fluidization has developed. Maximum 
effect was observed at sonic re- 
sonance for low frequency sound 
(200-600 cycles per second). 

The work described is exploratory; 
attention is now being directed at the 
attenuation in beds of solids, the 
suitable design of reactor cavities, the 
effect on surface reactions, and higher 
frequency effects. Whether the ex- 
perimental conditions described have 
counterparts in practice has not yet 
been ascertained. 


Sphere dia. = 0.313 in, 
Air flow 1.0 c.f.m. 


Duration 300 min. 


---- Quiet run 


mee 376 6.9:0., 160 db. 


Figure 8—Typical profile of single 
sphere of naphthalene showing the 
effect of sound upon evaporation. 


Nomenclature 


A = sphere surface area, sq. ft. 
1), = sphere diameter, ft. 
D, = diffusivity of naphthalene vapour 
in air, ft.? hr.-! 
= voidage of bed 
G = mass velocity, superficial, lb. hr. 
it. 
= mass velocity, superficial molar, 
lb. mole hr.~! ft.-? 
mass transfer factor 


1 


Kg . pam . (Ns-)?/3 


S/A . in (1 — poe/ps)~! (Nse)?!8 


mass transfer coefficient, lb. moles 
hr.~! ft.-? atm.-! 
AW Mast st. PLM 
molecular wt. of naphthalene, Ib. 
Ib. mole 
air viscosity, lbs. hr.~! ft.-! 

, D,.G 
particle Reynold’s Number, 


modified Reynold’s Number, 
D,.G 


u(1l — €) 


Schmidt number, . 

pDy 
mean partial pressure of air in 
interfacial film, atm. 
log mean partial pressure of naph- 
thalene in air, atm. 

p air density, Ib. ft.~* 

S tube cross section, ft.? 

t duration of transfer period, hrs. 

Vs superficial air velocity, ft. hr.~! 

AW wt. loss of naphthalene from the 
column, lbs. 
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Studies on Flocculants for 


Settling, Thickening and Filtration 


in the Sherritt Gordon Process’ 


D. J. I. EVANS?, 


W. KUNDA*, 


Excessively slimy leach residues were encountered during the first 


few months of operation of the Sherritt Gordon Process. 


The effect of 


these slimes on the liquid-solids separation of the leach liquor from the 
residues is described and the laboratory investigations to find the most 
suitable flocculant for improving the characteristics of the slurries are 


detailed. 


The paper is concluded with a discussion of the effect of the 


flocculant selected, Separan 2610, on the refinery operation. 


HERRITT Gorpon Mines LIMITED 
S employs an ammonia pressure 
leaching process in their refinery ‘!) 
at Fort Saskatchewan, Alberta, for the 
recovery of nickel, copper and cobalt 
from the nickel concentrate produced 
at Lynn Lake, Manitoba. 

The process (2: 3: +. 5) involves leach- 
ing the sulphide concentrate under 
pressure with air and ammonia to 
dissolve the nickel, copper, cobalt and 
most of the sulphur, and to leave the 
iron and siliceous material as a dis- 
cardable residue. The copper is re- 
covered from the solution as a high- 
grade sulphide precipitated during 
distillation stage for the re-circulation 
of part of the ammonia. After further 
solution purification, to remove traces 
of copper and undesirable sulphur 
compounds, the nickel is precipitated 
by hydrogen under pressure. The 
cobalt, together with unreduced 
nickel, is precipitated as a mixed 
sulphide (by hydrogen sulphide) for 
eventual cobalt recovery and the end 
solution is ev aporated to produce 
crystalline ammonium sulphate. 

The leaching operation (the circuit 
is shown in Figure 1) is essentially a 
continuous, two-stage, counter-cur- 
rent system with fresh concentrate 
being partially leached (60%) in the 


leach liquor produced in the second 
stage leaching circuit (adjustment and 
final leach — ely). Such a sys- 
tem is necessary o produce a solu- 
tion suitable for. treatment in the 
copper sulphide precipitation step. 

Horizontal, water-jacketed auto- 
claves, fitted with four agitators 
operating in compartments separated 
by overflow weirs to give a certain 
amount of staging during leaching, 
are used in both leach circuits. The 
adjustment leach (two autoclaves in 
parallel operating at 170-180°F and 
90-100 psig.) is followed by a thick- 
ener-filter, liquid-solids — separation 
circuit. The final leach, two parallel 
trains of three autoclaves processing 
(under the same conditions) the re- 
pulped partially leached solids, is also 
followed by a thickener-filter liquid- 
solids separation circuit. In the latter, 
three washing stages, consisting of 
repulping and refiltering, are given to 
the residues to ensure a maximum 
recovery of soluble metal values be- 
fore discarding the leach residues to 
a pond. 

The thickeners are 105 ft. diameter, 
completely enclosed. Hardinges 
thickeners; the filters are Dorr-Oliver 
American Disc filters with seven discs 
on the adjustment filter and ten discs 


NOTE: This article is based on a paper presented at the Third Western Regional Conference of 
The Chemical Institute of Canada, ae September 13-15, 1956. 


1Manuscript received February 17, 
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2 Assistant Head, Research Department, Sherritt Gordon. 
%Senior Research Chemist, Research and Development Laboratories, Sherritt Gordon. 
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on the four filters in the final leach 
circuit. 

A typical analysis of the nickel 
concentrate (produced by ae 
flotation) is given in Table 1. The 
mineral species identifiable in the 
concentrate under the microscope 
are: 

Pentlandite, (Fe,Ni)S 35-49% 

Chalcopyrite, CuFeS,, 3-6% 

Pyrite, FeS, 1-4% 


Py rrhotite, Fe. vs 24-42% 
®Where x varies be tween 0 and 0.2. 


together with various siliceous ma- 
terials such as tale and quartz. 
During leaching it is assumed ‘) 
that the sulphur in the pentlandite, 
chalcopyrite and part of the pyrrho- 
tite is successively oxidized through 
the thiosulphate and _polythionate 
species to sulphate and sulphamate. 
Simultaneously the nickel, copper 
and cobalt are converted to soluble 
ammine complexes and the iron in the 
pentlandite, chalcopyrite and part of 
the pyrrhotite is oxidized to a spongy 
hydrated iron oxide. The latter is 
readily broken down by the agitation 
necessary for satisfactory leaching. 
During the first few months of 
operation, difficulties were ex- 
perienced in the liquid-solids separa- 
tion circuits due to an excessive slimes 
load. Briefly the experience was as 
follows: satisfactory continuous oper- 
ation would be achieved for a few 
days although during this period slow 
settling pulp would be discharged into 
the final thickener: the slimes level 
in this thickener would gradually rise 
and eventually overflow with the 
solution recycled to the adjustment 
leach circuit: due to the recycled 
slimes, extremely slow settling pulp 
would then be discharged to the ad- 
justment thickener and the slimes 
level in this thickener would rise to 
overflowing: the remainder of the 
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Final Leach Solution 


Nickel Concentrate 
From Rail Cars 


Screen and Crusher 


Repulp Tank 
Spent Air to 
H 175°F Ammonia Recovery 


iy Adjustment Leach 
100 psi Air Circuit. 


Autoclaves 


Slurry Coolers 
and a 
Pressure Release 5 


Valves 7 


Adjustment 
Thickener Adjustment Leach Solution—3 


to Copper Boil, etc. 
Adjustment Filter 
epulp Tank *“——————Aqueous Ammonia Recycle 


(from Ammonia Recovery 
Circuit ) 
Ammonia 
Air 


Final 
Leach 


Auto- 175°F 


Cleves C+ 100 psi Air 


Slurry Coolers 
and 
Pressure Release 
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Final Filter 

Repulp Tank Aqueous Ammonia Recycle 
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Repulp Tank 
2nd Wash Filter 


Repulp Tank 


3rd Wash Filter 


Repu [p Tank Water 


Leach Residues to Pond 


Figure 1—Leach circuit with liquid-solids separation. 
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circuit, not being designed to cope 
with large quantities of semi-colloidial 
hydrated iron oxide, would cease to 
function satisfactorily and the con- 
ceutrate feed rate to the circuit would 
have to be cut to a minimum until the 
conditions in the thickeners improved. 

An intensive study of the plant 
operation, with a correlation of the 
settling characteristics of the leach 
slurries in various parts of the circuit 
with the operating variables and pro- 
cedures, was undertaken. This study 
showed that the pulps were charac- 
terised by progressively poorer set- 
tling rates as the solids passed through 
the circuit, indicating that the leach 
process had an important influence 
on the settling characteristics. It also 
showed that extremely slow-settling 
slurries were produc ed when the feed 
to the autoclaves had to be shut off 
and agitation (to prevent sanding up 
of the autoclaves) was maintained 
without leaching or discharge of the 
solids. Operation at low feed rates 
was also found to produce slow- 
settling slurries. Thus the indications 
from the plant study were that over- 
agitation was the main factor contri- 
buting to the extremely slimy nature 
of the leach residues; this was proved 
conclusively later and steps were 
taken in the plant to provide adequate 
agitation without excessive degrada- 
tion of the leach residues. 

One aspect of the plant study in- 
cluded the re-determination of the 
thickener area requirements ‘®7) 
since the original design data was 
obtained from the Ottawa pilot plant 
where no difficulties had been ex- 
perienced in the liquid-solids separa- 
tion. It was shown that the thickeners 
available would have sufficient ca- 
pacity (without allowing the standard 
30% safety factor in design) to 
handle 200 tons per day of concen- 
trate under normal conditions but 
with slow-settling pulps (such as those 
discharged after a temporary shut- 
down for maintenance and those pro- 
duced when the slimes were re-cycled 
to the head of the leaching circuit) 
the capacity would be reduced to 
below 100 tons per day of concen- 
trate. The desired throughput, at that 
time, was 200 tons per day of con- 
centrate. 

Thus the necessity of improving 
the settling and thickening character- 
istics of the leach slurries was evident 
if considerable capital expenditure on 
further thickeners and filters was to 
be avoided. The use of flocculants 
sugested itself as an obvious and im- 
mediate remedy. 


Flocculant behaviour 
The chemical, screen and sub-sieve 


TABLE 1 


NICKEL CONCENTRATE, 
CHEMICAL AND SIZE 


Chemical Analysis 
INE. 
G6: 
Cu. 
So 
Insol. 
Ce. 


Calculated Mineralogical Composition 
Pentlandite, (Fe,Ni)S 
Chalcopyrite, CuFeS, 

Pyrite, FeS» 

Pyrrhotite, Fe;_xS. 

Ferric Oxide Fe.O;. 

Adsorbed Ni (on FesO3nH:;O).. 


(Assumed, 


Nickel 
Concentrates 


53% 


ADJUSTMENT AND FINAL LEACH RESIDUES 
ANALYSES, Nov.-Dec. 1954 


Final 
Leach 
Residues 


Adjustment 
Leach 
Residues 


34% .28 


0.42 “42 :29 


49 24 13 
. 86 40.70 7.80 
.86 22.23 wae 
.78 5.80 7.30 


0.19 18 


4 


34. 


nil (assumed) 


2.50% 13% 
4. 21 
2.02 2.4 


31 0.16 
(assumed) (assumed) 

15 30.56 33.70 

16.50 37.60 

not determined 0.64 


Fe/ Ni ratio in Pentlandite = 1.0, 


Assumed x in Pyrrhotite formula =O, 
Degree of hydration of Fe.O; not known). 


Screen Analysis (Tyler) 


+ 65 mesh. 
~65 + 100 
—100 + 150.. 
—150 + 200.. 
—200 + 325.. 
—325 


Sub-Sieve Size 
—40 + 30 microns 
—30 + 15.. 
—-15+ 5.. 


—§ 


Average S.G. of Solids 


size analyses (as determined by the 
beaker decantation method (®) of the 
solids to be separated in the two 
liquid-solids separation stages are 
given in Table 1. The slurries are 
20-25% solids by w eight. The liquors 
(s.g. 1.1-1.2), in both cases, are highly 
ammoniacal (pH 8-10) and contain 
ammonium sulphate, nickel, cobalt 
and copper ammine complexes and 
small quantities of sulphamate, poly- 
thionate and thiosulphate radicals, 
giving a total electrolyte concentra- 
tion of 450-500 g.p.l. 

The problem presented was to find 
the most economic addition agent 
that would improve the thicken- 
ing characteristics of these slurries 
adequately without affecting the 
chemistry or operation of the remain- 
ing stages in the Sherritt Gordon 
process. 

Although little is known concern- 
ing the mechanism of dispersion and 
flocculation it is possible to formulate 
a hypothesis based on the fundamen- 
tal physico-chemical principles gov- 
erning the behaviour of solid suspen- 
sions in aqueous solutions to account 
for the majority of the observations 
on the ammonia leach slurries. 
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“© weight of total solids). 


34. 32! 

9. 
£2. 10. 
12. 19. 
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All minerals, when suspended in an 
aqueous solution, are soluble to some 
extent and it is the tendency of the 
ions on the surface of the particles 
and within the liquid boundary layer 
to establish a condition of equilibrium 
that results in an electrical charge 
(zeta-potential) on the particles. The 
generally accepted concept of the 
electrical charge is of two diffuse 
layers of oppositely charged ions (the 
Helmholtz and Gouy ‘layers). Al- 
though it has not been proven 
(streaming potential studies would be 
required to verify this point) it is 
probable that the hi igh electrolyte 
concentration of the ammonia leach 
slurries, with the correspondingly 
marked tendency for ionization, re- 
sults in a high electrical charge on the 
suspended particles. The metal ions in 
the sulphide minerals, the absorbed 
nickel on the hydrated iron oxides (19) 
and the metal ions in solution are the 
ions common to both solid and liquid 
phases. 

It has been shown experimentally 
(11) that it is the zeta-potential on 
the suspended particles that primarily 
controls the slime-coating on mineral 
particles and it is suggested that the 
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Figure 2—Settling rate curves, with flocculants. 


A—Final leach discharge 


B—Final leach discharge with 2.0 lbs./ton*, causticized Beetsol starch 
C—Final leach discharge with 2.0 lbs./ton*, uncausticized Beetsol starch 


D—Final leach discharge with 0.62 lbs./ton*, 
E—Final leach discharge with 0.06 Ibs./ton*, 


Aerofloe 552 
Separan 2610 


F—Final leach discharge with 0.31 lbs./ton*, Separan 2610 


“Of Solids. 


conclusions drawn from this work 

can be applied equally to the mechan- 
ism of flocculation. Thus flocculation 
is governed by: 


1. The sign and magnitude of the 
charge on the suspended particles. 
Thus flocculation is facilitated by 
the mutual attraction of oppositely 
charged particles in a multi-solid 
phase system and is inhibited by 
the electrostatic repulsion of simi- 
larly charged particles. 

2. The surface tension between the 
liquid and solid phases. This is an 
influencing factor only when the 
particles have actually made con- 
tact with one another. 

3. The Van der Waals forces govern- 
ing the attraction of solid particles. 
Here again these forces are weak 
and are only of significance when 
the interpenetration of the outer 
ionic layers has been achieved. 


4. The average particle size and the 
particle size range ‘'?); the form 
and shape of the particles; and the 
form, size and shape of the result- 
ing aggregates. These factors are 
also of less significance than the 
electrical charge on the particles. 
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The solids concentration in the 

suspension. 

With regard to the electrical charge 
there is a certain critical zeta- poten- 
tial above which dispersion is main- 
tained and below which orthokinetic 
coagulation and _ flocculation — takes 
place (13), This critical zeta-potential 
is dependent upon the physical and 
chemical nature of the solid particles, 
the solids concentration and the ionic 
atmosphere of the liquid phase. Be- 
low the critical charge the particles 

can approach one another sufficiently 

for the surface tension and Van der 
Waals forces to become effective and 
cause flocculation ('), With only the 
latter forces binding the particles to- 
gether in the floccules, the bonds may 
be readily broken by mechanical 
forces. 


Laboratory testing of flocculants 

Confirmation that the solids in the 
ammonia leach discharge _ slurries 
were well dispersed was afforded by 
the fact that the sands settled rapidly 
away from the slimes (i.e. free set- 
tling conditions were prevailing). 
With a flocculated pulp the heavier 
and coarser particles would have been 
bonded with the slimes and the re- 


sulting aggregates, 
specific gravity unleached sulphides, 


due to the high 


would have settled much more 
rapidly than the slimes alone. 

Thus to prevent the inherent dis- 
persion of the slurries (and since it 
was impossible to reduce either the 
(OH)>~ ion or electrolyte concentra- 
tions) it was necessary to find a ma- 
terial that would provide either a 
coagulating ion (such as Cat+++) or 
a protective coating on the particles, 
preventing their ionization tendency 
and reducing the zeta-potential to be- 
low the critical. 

No effective coagulating ion could 
be found; lime, calcium salts, phos- 
phates and polyphosphates were tried 
without success. 

A considerable amount of work 
has been carried out to improve the 
sedimentation of the red muds in the 
Bayer Process (a similar system to 
the one discussed here) with the con- 
clusion that organic materials such as 
flour and starch were the most useful 
for improving the liquid-solids separ- 
ation of the colloidal iron oxide slimes 
from the caustic aluminate liquors ‘!5), 
Recent work has confirmed that the 
flocculating effect of these materials 
is accompanied by an almost com- 
plete neutralization of the zeta- -poten- 
tial on the particles’), For this 
reason the initial laboratory tests on 
the ammonia leach slurries were car- 
ried out with a wide range of starches 
and starch products. A number of 
commercial organic flocculants, avail- 
able at the time of the investigations, 
were included in the testing program; 
it was assumed that these flocculants 
would act like the starches. 

For the initial evaluation of these 
materials, sedimentation tests were 
carried out in one-litre graduated 
cvlinders on final leach discharge 
slurry samples treated with various 
quantities of the flocculants. On set- 
tling, the slurries gave a distinct de- 
marcation between the clear liquor 
and the settling solids. The depth of 
clear liquor was determined at regu- 
lar intervals; the settling rate (ac- 
tually of the slowest settling particles) 
was determined by plotting the depth 
of clear liquor to a time abscissa and 
calculating the gradient of the free- 
settling portion of the curve. Ty pical 
settling rate curves are given in 
Figure 2. Curves C, D and F show 
both feos settling (a to b), in which 
the particles fall unhindered by each 
other or at least where the elimina- 
tion of water is by a regular sub- 
sidence of the slimes level, and com- 
pression settling (c to d) in which 
the particles or flocs rest on, or are 
partly supported by, each other. Sec- 
tions b to c on the curves represent 
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the transition between free and 
compression settling. Curves A, B, 
and E are in the free-settling zone 
oniy; compression settling conditions 
would be reached if the data beyond 
the three hours were included. 

‘The settling rates of the final leach 
discharge samples taken for the 
flocculant testing series were far from 
constant as they were, necessarily, 
taken over a period of several weeks. 
For this reason the criterion used for 
evaluating the effectiveness of a 
flocculant was the settling rate ratio— 
the ratio of the settling rates, with 
and without the flocculant, on the 
same slurry sample. Good consistency 
in the results was obtained by using 
this criterion. 

Simultaneously with the sedimenta- 
tion tests, filtration tests were carried 
out using a 0.1 sq. ft. Eimco labora- 
tory test filter leaf covered with 
Duplan 2080 nylon filter cloth. The 
filter leaf was immersed in the 
lightly-agitated, conditioned slurries 
and the volume of filtrate collected in 
a given time was recorded. Here 
again the filtration rate ratios (with 
and without flocculants as before) 
were preferred to the actual filtration 
rates for comparative purposes. 

The settling rate, and filtration 
rate, ratios for some of those ma- 
terials which improved the aneeng 
characteristics are given in Table 2. 
The program included testing all the 
materials received on the final leach 
discharge slurries, selecting those 
which offered improvement and con- 
firming this tendency on the adjust- 
ment leach discharge slurries. Those 
giving the best results in both series 
were then further investigated. 

Of the commercial starches tested, 
Cellodex and W.W. Dextrine gave 
the best results when added in the 
causticised state at 2.0 lbs. per ton of 
suspended solids. Beetsol and Her- 
cules Regular starches, proved more 
effective when added in the uncausti- 
cized state and also more effective 
than the causticized Cellodex starch. 
The behavior of the starches was very 
specific and few of them brought 
about any improvement on the set- 
tling rate. 

Of the 30 commercial flocculants 
tested, Separan 2610 proved superior 
to all of the other flocculants, even at 
concentrations as low as 0.18 lbs. per 
ton of solids; it also improved the 
filtration rate more than any other 
addition agent tested. Reagent S-3000 
and Aerofloc 552 followed Separan 
in efficiency. 

The beneficial effect of Aerofloc 
552 had been shown early in the in- 
vestigations by the leach plant operat- 
ing staff and, being readily available, 
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it was added to the leach discharge 
slurries in quantities of 0.5-1.0 Ibs. 
per ton of concentrate feed, over a 
period of three months. Its improve- 
rate was ap- 
parent although sliming conditions in 
the thickeners still occurred at fre- 
quent intervals. With Separan 2610 


ment on the settling 


becoming available 
quantities a 
was recommended. 


Separan was 


in 


switch to this 


preferred 


commercial 
material 


because 


smaller quantities were required to 


give the 


same improvement in the 


settling characteristics so resulting in 
a lower reagent cost, less dilution of 
the leach liquors and a greater flexi- 


bility 


in application. 
also superior to the other 


Separan 


was 
flocculants 


in that the incidence of compacting 
settling took place at a higher pulp 
density (advantageous in a thickener 


employi ing a raking mechanism). 


compacting settling 


density 
could be obtained i 


SEDIMENTATION 


No. Addition 


1. | Blank 
2. | Cellodex Starch 


3.) W.W. 20 Dextrine 


4. Cato Starch 


Beetsol Starch 


wn 


6. Hercules Regular 
Starch 


Beetsol Starch 


~ 


8. | Hercules Regular 
Starch 


9. | Aerofloc 548 
10. | Aerofloc 552 
11. | S-3000 


12. | Separan 2610 


13. High Viscosity 


Cellulose Gum No. 


14. Hallmark, 
Jaguar MD-C 
15. | Acrysol GS 


16. | Lytron X886 
(also 887 and 889) 


The 


characteristics 
were also improved by Separan, in- 
dicating that the desired high pulp 


in the thickener 
a shorter 


TESTS ON FINAL 
Note: 


Re 


w 


20 


underflows 
time. 


TABLE 
LEACH DISCHARGE 


Flocculant Costs, 


It was also found that Separan and 
the starches, in addition to improving 
the filtration rate, gave an easily de- 
tachable filter cake whereas the un- 
conditioned slurry, and the Aerofloc 
552 and S-3000 conditioned slurries 
gave filter cakes which were firmly 
adherent to the filter cloth. 

Prior to recommending Separan for 
plant use, it was necessary to deter- 
mine its effect on the other process 
stages and to determine the most 
suitable point of application. It was 
apparent that the quantity to be 
added and the operational control 
could best be determined (using the 
laboratory figures as a guide) from 
observation of its effect on the plant 
oper ation. 

In the initial laboratory investiga- 
tions on the effect of Separan, an 
excess (5 Ibs. per ton) added to 
leaching charges in the laboratory 
autoclave, was shown to have no 
effect on the metals extractions; the 
settling rates after leaching were also 
unaffected. 

An excess (1.2 lb. per ton) was also 


) 


With Appition AGENTS 


January, 1955. 


Starches 2—6 added in causticized state 


(0.2 Ibs. 


Supplier 


Ogilvie Flour 


Ogilvie Flour 


National Adhesives 


Canada Starch 


Canada Starch 


Canada Starch 


Canada Starch 


Cyanamid 
Cyanamid 
Cyanamid 


Dow Chemical 


Hercules 
Stein-Hall 
Rohm and Haas 


Monsanto 


Caustic Soda per Ib. of starch) 
. Starches 7 and 8 added uncausticized. 


Quantity Ratios 


Ibs/ton 
Cost of : 
$/Ib solids | Settling Filtration 
Rate Rate 
1.00 1.00 
0.115 1.0 3.51 
0.115 2.0 6.04 
1.0 Ruan 
2.0 5.88 
1.0 1.45 0.93 
2.0 2.09 1.02 
0.115 1.0 3.07 1.40 
0.115 2.0 cna 2.04 
0.115 1.0 3.37 4:25 
0.115 2.0 3.80 1.88 
0.115 1.0 3.65 1.46 
0.115 2.0 8.70 1.78 
0.115 1.0 3:45 1.40 
0.115 2.0 7.50 1.92 
0.70 0.31 auo 0.88 
0.70 0.62 4.8 1.03 
0.75 0.31 Aco 1.08 
0.75 0.62 sad 1.31 
1.00 0.31 5.0 1.44 
1.00 0.62 9.1 1.61 
1.00 0.06 3.2 1.29 
1.00 0.18 9.5 2.09 
1.00 0.3 26.0 2.20 
1.00 0.62 80.0 2.78 
0.3 3.4 1.18 
0.62 4.0 1.63 
0.50 1.0 1.08 
1.00 1.0 0.97 
0.50 Zen 1.03 
1.00 4.2 1.03 
0.50 54 1.80 
1.00 4.5 2.34 










(settling rate without flocculant ) 
(settling rate with flocculant ) 


Settling rate ratio 


0.20 0.40 


Floecculant, Ibs./ton of solids 


Figure 3—Comparative flocculating efficiencies on final leach discharge (initial 
settling rate without flocculant—0.49” /hr. ) 


A. 
B. 
C. 


0.60 


Separan 2160 
Reagent S-3000 
Aerofloe 552 


D. Beetsol and Hercules regular starches (uncausticized ) 


TABLE 3 





RELATIONSHIP BETWEEN Dry SOLIDS FILTERED PER HOUR AND THE SEPARAN 
ADDITION TO A FINAL 


(P.D.—1.840) 
A. LABORATORY TESTS 
Separan 
Test No. Ibs./ton © Moisture 

of solids in Wet Cake 
a nil 19.5 1 
zs 0.05 21.5 1 
a 0.10 25.3 16. 
4. 0.15 23:2 20. 
>. 0.24 22.0 21 
6. 0.32 21.2 28. 
7. 0.40 20.3 | 31. 


B. Pitot PLANT FILTER TEsTs: 


On Slurry Already Thickened with 
1. nil 
2. 0.05 
3: | 0.06 
4. 0.10 
5. | 0.12 
6. 0.15 
7. | 0.24 
&. 6.25 
2. 0.48 

10. 0.96 


0.06 


19. 
21. 
22 
22 
23 
23. 
22 


SANOHK OHNO 


Ibs.of dry 
solids per sq. 
ft. per hour 


4. 
4. 


On WON 


THICKENER UNDERFLOW 


_ 


°> Improvement 
(based on Filtration 
Rate ratio) 


15.0 
41.9 
50.2 
95.5 
114.0 


Ibs. of Separan per ton of Solids. 


00 00 CO = = 00 00 


nA 


wo 


54 


18 


.90 
.00 
.69 


90 


.45 
“ol 
.03 
.00 


25 
36 
68 
78 
36 
29 
27 
38 
deteriorated by 
54% 


added in turn to the adjustment leach 
discharge slurry and solution, the 
copper-free solution, and the nickel 
reduction feed solution. The solutions 
were then processed through the re- 
maining stages. From this it was con- 
cluded that Separan added to the cir- 
cuit before the oxidation and hydro- 
lysis stage '-2) had no detrimental 
effect on the process. Separan added 
after that stage (at least in this 
quantity) caused flocculation of the 
nickel with consequently poor opera- 
tion in the’ reduction  densifica- 
tions (17), 

It was also shown by a solution 
switching procedure and correlation 
with the settling rates, that 90-100% 
of the Separan added to a leach slurry 
would pass out of the circuit with the 
solids. 

These experiments cleared Separan 
for use in the refinery circuit. 

To ascertain the most suitable 
point of application for the flocculant, 
series of stability tests were carried 
out on the six most promising floccu- 
lants, (Separay, S-3000, Aerofloc 552, 
Beetsol starch, Hercules Regular 
starch and Cellodex starch). It was 
shown that agitation alone readily 
caused breakdown of the flocculating 
bonds of Separan, S-3000 and Aero- 
floc 552, that agitation at elevated 
temperatures (175°F) accelerated this 
breakdown and that neither agitation 
nor temperature had any effect on 
the flocculating effect of the starches. 
From this work it was recommended 
that if Separan were used it should 
be added to the cooled discharge 
slurries at a point where the agitation 
was sufficient only for mixing prior 
to thickening. A two-point addition 
was adopted for each thickener feed. 


Effect of Separan on the process 
Thickening—Opportunity was 
taken to confirm the laboratory re- 
sults and conclusions during the first 
few months of adding Separan to the 
liquids-solids separation circuits. 


The relationship between the quan- 
tity of Separan added and the im- 
provement to the final leach residues 
settling rate is shown in Figure 3. 
Comparative relationships for S-3000, 
Aerofloc 552 and the best of the 
starches are also included. The advan- 
tages of Separan over the other 
flocculants are obvious. 


In earlier investigations it had been 
shown that settling tests at various 
pulp densities were necessary to pre- 
dict the thickener diameter require- 
ments. Thus it was necessary to con- 
firm that the faster initial settling 
rates with Separan would also mean 
improved thickening characteristics. 
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Figure 4—Thickener diameter requirements against Separan addition to the final leach discharge to give a final thickener 


Thickening tests ‘* *) with 0.0, 0.03, 
0.06, 0.10, 0.15 and 0.20 lb. of Separan 
per ton of solids were carried out on 
a sample of final leach discharge of 
good settling characteristics (1.0 in. 
per hour). 

The thickener diameters required 
for various feed tonnages are plotted 
to a Separan addition abscissa in 
Figure 4. The relationship between 
the initial settling rate and _ the 
Separan addition is included in the 
inset. The curves show that, without 
Separan, a concentrate feed rate of 
200 tons per day would require a 
103 ft. dia. thickener in the final leach 
circuit; this would leave a safety fac- 
tor of only 1.9% for abnormal 
operating conditions. With 0.03 Ibs. 
of Separan per ton of concentrate the 
diameter required would be reduced 
to 73 ft. giving a safety factor of 
43.7%. This is a more than adequate 
safety factor. Alternatively, the ad- 
dition of 0.03 lbs. of Separan per ton 
would allow the residues from a con- 
centrate feed of 250 t.p.d. to be 
separated in the 105 ft. thickener 
available and still leave the 30% 


underflow P.D.—1.850. 


safety factor specified normally in 
design calculations. 

Increasing the Separan above 0.03 
lb. per ton progressively improves 
the thickening characteristics of the 
slurry, although not by as much as 
does" the initial addition. Thus the 
optimum addition for up to 300 tons 
per day rate would aparently lie be- 
tween 0.03 and 0.06 lb. of Separan per 
ton and the inferences drawn from 
the settling rate data would appear to 
be justified. 


Filtration—In addition to improv- 
ing the thickening characteristics 
Separan also improved the filtration 
characteristics of the thickener under- 
flows. A laboratory series to illustrate 
this, using the Eimco Corporation 
test procedure ‘!8) is given in Table 
3. Final thickener underflow was used 
for the filtration tests and Duplan 
2080 (the standard cloth for the plant 
filters at the time) was used on the 
laboratory test filter leaf. In Figure 5 
a comparison of the effect of Separan 
on an underflow sample taken under 
normal operating conditions (A) is 
made with one taken when an exces- 
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sively high slimes burden was present 
in the leaching circuit (B). 

The curves show that above 0.05 
lb. per ton there is an approximately 
linear relationship between the 
Separan addition and the filtration 
rate (as determined by the weight of 
dry solids per sq. ft. per hour). The 
fact that more Separan is required on 
a slimy pulp to achieve the same 
filtration rate is also illustrated. From 
the figures reported, Separan would 
appear to increase the moisture con- 
tent of the filter cake at low con- 
centrations. In later studies this trend, 
although definite, was shown to have 
little significance. 

Pilot plant filter tests, using a two 
compartment, two-disc Eimco filter, 
were carried out to confirm the 
laboratory results. Using this filter, a 
side-by-side comparison on the same 
final thickener underflow slurry, with 
and without a flocculant, was pos- 
sible. The filter media on both discs 
was Duplan 2080 (a multifilament 
nylon cloth). 

The results from these series (one 
series is reported in Table 3) showed 
a definite optimum Separan addition 
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Filtration Rate, Ibs. dry solids/sq. ft./hr. 


0.1 0.2 


0.3 0.4 


Separan addition, lbs./ton of solids. 


Figure 5—Laboratory tests, effect of separan on filtration rate of final thickener 


underflow. 
A. Normal slurry 


B. Slurry from thickener during sliming conditions. 
C. Filtration rate needed to process residues from 250 tons of concen- 
trate/day on existing 8'6” x 10 dise filter (assuming 30% safety factor 


and 100% cake discharge). 


for improving the filtration of 0.10- 
0.12 Ib. per ton (in addition to the 
0.06 lb. per ton previously added to 
the thickening stage). At this con- 
centration an improvement of 60-80% 
was obtained. Thus there would ap- 
pear to be a contradiction between 
the laboratory results, giving a linear 
relationship, and the pilot plant filter 
tests, giving an optimum, for the 
filtration efficiency versus Separan 
addition. This anomaly may be ex- 
plained by the more efficient agitation 
in the laboratory tests giving a more 
uniform cake build-up and the 100% 
cake discharge in the laboratory tests 
compared with 30-40% from the disc 
filters under the conditions of the 
test. With the higher Separan addi- 
tions, consolidation settling of the 
slurries in the filter boots took place 
and this was a further factor pre- 
venting the pick up of a uniform cake. 

The pilot plant results were later 
confirmed in the commercial plant. 
Since that time a combination of 
other improvements, to improve the 
uniformity of the cake build-up and 
the ease of cake discharge together 
with the addition of Separan, has im- 
proved the capacity of the tailings 
filter circuit well beyond the design 
capacity. Thus the need for extra 
filters, as the production rate of the 
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refinery was gradually increased, was 
avoided. 

Start-up operations for the refinery 
commenced at the end of May, 1954 
with each stage being brought into 
operation successively “) and with 
the first nickel powder being pro- 
duced on July 21. Between that date 
and the end of January 1955, when 
the full capacity of the refinery was 
attained, difficulties were frequently 
experienced in the leach plant due to 
sliming conditions throughout the 
leaching circuit. 

Separan was received in commer- 
cial quantities at the refinery on 
January 19, 1955. At that time opera- 
tion of the adjustment thickener was 
satisfactory with 30-40 in. of clear 
liquor, but sliming conditions existed 
in the final thickener with solids 
overflowing back to the adjustment 
leach. From the laboratory investiga- 
tions, Separan additions of 0.03 Ib. 
per ton of concentrate to the cooled 
adjustment leach discharge, 0.06 Ib. 
per ton to the cooled final leach dis- 
charge and 0.12 Ib. per ton to the 
final thickener underflow were re- 
commended. Two-point additions of 
a 0.25% Separan solution in cold 
water were employed. 

By January 24, 70 in. of clear liquor 
were being maintained in the adjust- 


ment thickener and by January 28 
the steady improvement in the depth 
of clear liquor (5 in. per day) had 
resulted in 45 in. of clear liquor in the 
final thickener. Apart from an initial 
over-loading of the rake mechanism 
when the slimes load was dropped 
too rapidly due to a considerable ex- 
cess of Separan, the operation of the 
thickeners has since been excellent. 

Thus from the observation of the 
plant operation, Separan (the effec- 
tive constituents of which are high 
molecular weight acrylamide poly- 
mer hydrolytes ‘!®)) has been shown 
to improve the thickening and filtra- 
tion characteristics of the ammonia 
leach residues and to have no harm- 
ful effect on the chemistry or opera- 
tion of the remaining process steps 
when added to the liquid-solids 
separation circuit in quantities of 
approximately 0.2 lb. of Separan per 
ton of nickel concentrate. 
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The Decomposition, Oxidation, Ignition 


and Detonation of Fuel Vapors and Gases 


Part 31. 


The Effect of Flow Configuration 


in an Annular Reactor of Vycor Glass 


on the Oxidation of Pentane 


R. O. KING?, 


The oxidation and ignition of n-pentane-air mixtures in an annulus 
1.0 mm. wide between concentric tubes of Vycor glass has been investi- 
gated in a flow system at atmospheric pressure and at temperatures 
rising to 800°C. The mixtures were supplied at rates ranging from 400 
to 50 cc./min. and the flow was necessarily laminar. Neither oxidation 
nor decomposition reactions occurred at surface temperatures lower than 
625°C. when a correct mixture was supplied at 400 cec./min. but when 
time of residence was doubled an oxidation reaction with a negative tem- 
perature coefficient occurred over the temperature range 470°C. to 
505°C. There was then a temperature range of 145°C. over which no 
reaction could be detected before the initiation of a high temperature 
raction. The maximum rate of the initial reaction increased and the 
temperature range of practically no reaction diminished as the rate of 
mixture supply was reduced to 100 cc./min. The products of the low 
temperature reaction were essentially aldehydes, acids and steam but 
there was no fluorescence, pressure fluctuation or ignition. It is shown 
that the negative temperature coefficient is a surface phenomenon. Thus 
it has not been necessary to attribute it in whole or in part to a gas phase 
reaction dependent upon a chain reaction mechanism. 


tors in subsequent studies of this 
phenomenon. Some claimed that the 
essential intermediate was a perox- 
ide ‘3) or a peroxide radical ‘*) while 
others attributed the results to the 


Introduction 


A negative temperature 
coefficient of the velocity of 
the low temperature oxidation of pro- 
pane and the butanes was first de- 


scribed by Pease who attributed it, at 
the time, to the inhibitory action of 
excess oxygen or diluent nitrogen at 
the higher temperatures ‘'). Neumann 
and Aivazov ‘2) used a bulb method 
to study the oxidation of n-pentane 
and attributed the negative tempera- 
ture coefficient of the primary re- 
action to autocatalysis by one of the 
intermediate products of which the 
maximum concentration was con- 
sidered to decrease with increase of 
temperature. This general explanation 
was elaborated by other investiga- 


accelerating effect on the initial re- 
action of various aldehydes, an effect 
which was found to diminish with 
increasing temperature in the low 
temperature range (5), (6), In any case, 
the general opinion was that the low 
temperature mechanism of the reac- 
tion was that of a degenerate explo- 
sion of the branched chain type and 
the rather sharp suppression of re- 
action with temperature increase was 
attributed to suppression of branch- 
ing by the further reaction of some 
important intermediate. 


!Manuscript received March 3, 1957. 
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S. SANDLER’, and R. STROM‘ 


A study of the references men- 
tioned and many others indicates a 
decided lack of agreement in respect 
of the mechanism of the chain reac- 
tion to which the negative tempera- 
ture coefficient was attributed. The 
magnitudes of the coefficient as 
shown by their experiments depended 
upon whether the flow or bulb 
method was used. Flow configuration 
conditions of the experiments dif- 
fered according to the dimensions 
and conformation of the reactors but 
its effect was not considered in dis- 
cussions of the cause of the negative 
temperature coefficient. 

Flow configuration is however 
shown by the experiments with pen- 
tane-air mixtures described in Part 2 
(7), to be a major factor in respect of 
reaction velocity and the magnitude 
of the negative temperature coeffi- 
cient. When, for example, pentane in 
a mixture with air was oxidized in a 
25 mm. diameter pyrex combustion 
tube, a negative temperature co- 
efficient was obtained when it was in 
the horizontal position but not when 
flow configuration was altered by 
using the tube in a vertical position, 
as shown by the graphs of Figure 1, 
reproduced from the reference 
quoted. 

Negative temperature coefficients 
were also obtained when King reac- 
tor No. 10 was used (8: P- 333), This re- 
actor was designed and used mainly 
to obtain conditions in which the 
oxidation of pentane went to final 
products when doped with iron 
carbonyl. The design of the reactor 
was such that oxidation of pentane, 
either with or without the dope, 
became a surface effect. Experiments 
were carried out with various rates 
of mixture supply and _ mixture 
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Figure 1—The change in the tempera- 
ture coefficient of reaction due to 
changing a tube reactor from the 
horizontal to the vertical position. 


strengths and are described in Parts 
1 (8), and 2 (7), and in Part 7 (©) from 
which results obtained when the dope 
was not used are shown by the repro- 
duced graphs of Figure 2. They show 
the increase in percentage of oxygen 
reacted in the low temperature re- 
gion that accompanied an increase in 
the concentration of pentane in the 
mixture with air. When the con- 
centration was double that required 
for a chemically correct mixture, the 
rate of oxygen reacted attained a 
maximum value of 33% at 415°C. 
The rate then diminished to 16% 
over the temperature range 420°C. 
to 500°C. The low and high tem- 
perature reactions were not separated 
by a temperature range of no reac- 
tion. Explosive ignition did not occur 
although the temperature was raised 
to 650°C. 


It was indicated by experiments 
with vertical combustion tubes de- 
scribed in Part 23 ‘!°) that aldehydes 
formed on the tube surface in the low 
temperature range were carried from 
the boundary layer into the central 
stream by circulating convection 
currents, and there accounted for the 
occurrence of fluorescence, pressure 
waves and the consequent mild type 
of non-explosive ignition. The oc- 
currence of ignition in the low tem- 
perature range thus depended upon 
the existence, apart from the boun- 
dary layer, of a central stream in 
which there were circulating con- 
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vection currents. Similarly, ignition 
in the high temperature range was 
attributed to the effect of the con- 
vection currents in removing oxida- 
tion products from the boundary 
layer in order to allow access of less 
reacted mixture to the surface. By 
analysis of the effluent from the tube 
no distinction could be made between 
the reactions that occurred on the 
surface and those that might have 
occurred in the central stream. It was 
decided, therefore, to eliminate the 
effect of the central stream by using 
a reactor in the form of a narrow 
annulus between two concentric 
tubes. The mixture flow through 
was necessarily laminar thus circulat- 
ing convection currents were elimi- 
nated and boundary layer conditions 
obtained. The results of experiments 
made accordingly are given in this 
Part. 


The term “negative temperature 
coefficient” is used to mean the dif- 
ferential quotient of the reaction 
velocity with respect to temperature, 
which can be positive or negative. 
The phenomenon is also but not com- 
monly, described as the factor by 
which the reaction velocity changes 
with an increase of 10°C. in temper- 
ature. By this arbitrary method, the 
factor can be greater or less than 
unity but is never greater than zero 
in conditions of experiment generally 
used. It cannot be applied to condi- 
tions used for the experiments de- 
scribed in this Part in which there 
is a temperature interval of zero re- 
action which is greater than 10°C. 
between the low and high tempera- 
ture reactions. 


Experimental 


The concentric tubes of the reactor 
were of Vycor glass. The inside 
diameter of the outer tube was 19.0 
mm. and the wall thickness 3.0 mm. 
The outside diameter of the inner 
tube was 17.0 mm. and the wall 
thickness 2.0 mm. The width of the 
annulus was thus 1.0 mm. Vycor 
glass is 96% silica. Tubes of pure 
silica were not available in the re- 
quired diameters. 

The design of the reactor and its 
arrangement in a vertical tubular 
electric furnace 2.0 in. in diameter 
are shown by Figure 3. Concentricity 
of the tubes was maintained by small 
bosses spaced 120° apart at the rela- 
tively cool ends of the reactor. The 
heated length of the furnace was 6.5 
in. and temperature as measured on 
the exterior surface of the outer tube 
attained a maximum value at approxi- 
mately the midsection of that length. 
This temperature is given in subse- 
quent text as related to rates of re- 


action and to the occurrence of 
ignition. There were steep tempera- 
ture gradients from the position of 
maximum temperature to the rela- 
tively cool ends of the reactor. On 
the other hand, the temperature 
variation across the reactor at the 
position of maximum temperature 
was always very small, being usually 
about 2 or 3°C. 


The methods used for the measure- 
ment of temperature, rates of air and 
pentane supplied, oxygen reacted and 
oxides of carbon formed, were as 
described in the appendix to Part 
23 (10, pp. 125-126), Unsaturates were 
determined, in one set of conditions 
only, by absorption in a saturated 
solution of mercuric sulphate in a 
25% sulphuric acid solution using the 
conventional Orsat apparatus. A 
modification of the methods used by 
Widmaier and Mauss “!!) in the 
analysis of the products of oxidation 
of gasolines was adopted for the de- 
termination of total aldehydes and 
peroxides. 


Chemically correct mixtures are 
described as “correct” or by the 
abbreviation C.C.M. and contain 2.63 
cc. of the vapor of pentane in 100 cc. 
of air. Mixtures are described as per- 
cent weak or rich according to the 


ain SUPPLY 200cz7h0] 


PENTANE SUPPLY, —| |/| 
A,100% RICH 


B, 50% RICH 
C, 25% WEAK aa 


aa 


PERCENTAGE OF AVAILABLE O, REACTED 


300 400 500 600 700 


TEMPERATURE °C. 


Figure 2—Oxidation of pentane in 

King reactor No. 10 and the effect of 

mixture strength on the _ negative 
temperature coefficient. 
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percentage by which the pentane 
vapor is less or more than that re- 
quired for a correct mixture. Tech- 
nical grade n-pentane as supplied by 
the Phillips Petroleum Company was 
used for the experiments. 


Flow Configuration in the 
Annular Reactor 


The Reynolds number for the flow 
of air through an annulus is calcu- 
lated on the basis of the mean hy- 
draulic radius being + of the dif- 
ference in the diameters of the 
boundary surfaces. Reynolds numbers 
thus calculated for the annular reac- 
tor, Figure 3, for a rate of air supply 
of 50 cc./min. were 1.91 at 25°C., 
1.19 at 300°C., 1.07 at 400°C. and 
0.98 at 500°C. 


The Reynolds number increases in 
direct proportion with increase in the 
rate of air supply. The maximum rate 
used for the experiments was 400 
cc./min, and they were always car- 
ried out at annulus temperatures ex- 
ceeding 300°C. The Reynolds num- 
ber was therefore never greater than 
9.60. This is so far below the critical 

value of approximately 2300 required 

for the initiation of turbulence that 
the effect of small percentage addi- 
tions of pentane vapor can be ne- 
glected and flow through the annulus 
1.0 mm. wide regarded as laminar in 
all of the experimental conditions. 
Thus there would always be a sta- 
tionery layer of gas on the boundary 
surfaces and the velocity profile of 
the flow would be parabolic. A dif- 
ference in the temperature of the 
boundary surfaces would cause 
merely a distortion of the parabolic 
profile, without creation of circulat- 
ing convection currents. 


Results 


1. The Effect of Residence Time, 
Correct Mixtures with Air 
Supplied at Rates of from 
400-50 ce./min. to the Annu- 
lar Reactor 

The continuance of reaction that 

begins on the heated surface must de- 
pend upon diffusion of the products 
from it in order to permit access of 
reacting mixture to it. The rate of 
surface reaction will therefore de- 
pend upon the rate at which the 
products escape from the surface and 
diffuse into the mixture _ passing 
through the annulus. The rate will 
thus depend upon the time of resi- 
dence of the mixture in the reactor. 
This time is basically proportional 
inversely to the rate of supply of the 
mixture which always contains an 
excess of diluent consisting of nitro- 
gen and unreacted oxygen. The range 
of residence times investigated is 


The Canadian Journal of Chemical Engineering, June, 1957 


AL.- CR. 


N BS 






WLM 


©0000 8800080800008 


ECLEPZLCEPLELELEEER 





GLEALEELELALEEALL EA 


POM 


+CR. 


7 ~~ c0tissvvsssicscsssisscdcccdkbbbbh 
6-5 ins 





a 
°o 
Oo 
oO 
Zz 
be 
=a 
= 
e LOWER END 

) d 

} 

° H y 

] y 

® y yi 

H y 

| y 

y 

RR) BH 

AI 1 

| 1 BY 

sd] 1K 

Ka x 


Figure 3—Reactor of Vycor annulus 1.0 mm. wide. 


from 0.4 to 3.2 sec. as calculated for 
an average gas temperature of 500°C. 
and corresponding to rates of air 
supply varying from 400 cc./min. to 
50 cc./min. 

Experimental results obtained to 
show the influence of temperature on 
rates of oxygen consumption and 
formation of oxides of carbon are 
given by the graphs of Figures 4 to 
8 for the several rates of mixture 
supply. 


(a) Air Supplied at 400 ce./min. 


At the shortest residence time used, 
Figure 4, measurable reaction did not 
occur until the maximum tempera- 
ture exceeded 625°C. The product 
of the reaction then contained the 
oxides of carbon and a small propor- 
tion of unsaturates as would be ex- 
pected if decomposition as well as 
oxidation of the pentane had oc- 
curred. The rate of formation of 
both carbon oxides increased rapidly 
with temperature. However, carbon 
monoxide production began at a 
lower temperature than that of car- 
bon dioxide and continued at a 








greater rate as the temperature was 
raised until the mixture ignited at 
745°C. when little more than 50% of 
the available oxygen was reacted. 
The yield of unsaturates reached a 
maximum some 20°C. before ignition 
occurred and then diminished rapidly. 


(b) Air Supplied at 200 ce./min. 


When the residence time was thus 
doubled reaction began at 420°C. as 
shown by Figure 5, and the rate in- 
creased to a maximum of 0.15 moles 
of oxygen reacted per mole of 
n-pentane, at 370°C. The temperature 
coefficient then became negative and 
reaction ceased to be measurable at 
505°C. and so remained at tempera- 
tures rising to 650°C. The high 
temperature reaction then began and 
proceeded with increasing velocity 
until ignition occurred at 745°C. 
when approximately 50% of the 
available oxygen was reacted. The 
oxidation products of the low tem- 
perature reaction were aldehydes, 
acids and water; carbon monoxide 
and dioxide formation being negli- 
gible. Those of the high temperature 
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annular reactor with pentane for a 
correct mixture with air supplied at 
400 cc./min. 
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Figure 5—Rate of reaction in the 


annular reactor with pentane for a 
correct mixture with air supplied at 


200 cc./min. 


reaction were the oxides of carbon 
and water. Unsaturates were also de- 
tected. Neither oxidation nor de- 
composition products could be de- 
tected in the intermediate tempera- 
ture range. A low-temperature frac- 
tional distillation (Podbielniak-Semi- 
Robot) of the effluent obtained at one 
point in this range indicated that 
n-pentane was the sole hydrocarbon 
present. 


(c) Air Supplied at 150 and 100 


ec. /min. 


Experimental results obtained when 
the rate of air supply was 150 
ce./min. are given by the graphs of 
Figure 6. Reaction began at 405°C. 
and proceeded with increasing ve- 
locity to a maximum rate of 0.70 
moles of oxygen reacted at a tempera- 
ture of 435°C. The temperature co- 
efficient then became negative and 
rate of reaction ceased to be measure- 
able at 585°C. The high temperature 
reaction began at 620°C. that is, after 
a temperature interval of 35°C. only 
of no measurable rate of reaction. 
Ignition occurred at 745°C., again 
when approximately 50% of the 
oxygen was reacted. 
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The experimental results obtained 
when the rate of air supply was 
further reduced to 100 cc./min. are 
given by the graphs of Figure 7. 
They differ from those obtained with 
a rate of supply of 150 cc./min. in 
that maximum rate of oxygen reacted 
per mole of pentane in the low 
temperature range increased from 
0.70 moles to 0.80 and remained 
nearly constant until the temperature 
coefficient became negative at 500°C. 
and reaction rate decreased to 
practically nil at 610°C. The high 
temperature reaction began 25°C. 
later, at 635°C. and as before pro- 
ceeded with increasing velocity until 
ignition occurred at 747°C. 


(d) Air Supplied at 50 ce./min. 

The experimental results are given 
by the graphs of Figure 8. A com- 
parison of the results w ith those pre- 
sented for the 100 cc./min. rate of 
supply is of interest. The following 
differences should be noted: (1) at 
the lower flow rate, reaction began 
at a temperature 50°C. lower, coin- 
cident with the formation of some 
carbon dioxide, and was followed 
about 15°C. later with the formation 
of the monoxide in addition to alde- 
hydes, acids and water as found i 
the products obtained at the lower 
residence time, (2) maximum rate of 
reaction in the low temperature range 
increased from 0.80 moles to 1.30 
moles of oxygen reacted, (3) the low 
and high temperature ranges of reac- 
tion were not separated by a tempera- 
ture range of no reaction thus 0.75 
moles of oxygen was reacted per 
mole of pentane at 625°C. whereas 
there was no reaction at that temper- 
ature when the air was supplied at 
100 cc./min. 


(e) Ignition Temperatures 

The ignition of a correct mixture 
occurred within the temperature 
range 745° to 750°C. with rates of 
supply ranging from 100 cc. to 400 
cc./min. Ignition was uncertain at 
rates less than 100 cc./min. and at 50 
cc./min. temperature was raised to 
800°C. without the occurrence of 
ignition, although the mixture 
strength was varied from 50% weak 
to 85% rich. 


(f) Production of Aldehydes 

The influence of temperature on 
the total aldehyde production is 
shown by Figure 9. A maximum was 
obtained at a temperature of 430°C., 
which corresponds to that at which 
the low temperature rate of oxygen 
consumption was a maximum. The 
rate of aldehyde formation then de- 
creased rapidly with temperature 
over the next 100°C. As oxygen con- 
sumption remained relatively con- 
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Figure 6—Rate of reaction in the 
annular reactor with pentane for a 
correct mixture with air supplied at 


150 ce./min. 
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Figure 7—Rate of reaction in_ the 
annular reactor with pentane for a 
correct mixture with air supplied at 


100 ec./min. 


stant over this range, see Figure 6, 
the formation of oxygen containing 
products other than those for which 
analyses were performed is thus in- 
dicated. 

The graph of Figure 10 shows the 
effect of increase of residence time 
on the rate of aldehyde formation 
when a correct mixture was used and 
the maximum temperature of the 
reactor was maintained at 430°C. It 
will be seen that the volume per- 
centage of aldehyde formation in- 
creased slowly as rate of air supply 
was increased from 50 cc. to 100 
ce./min. The maximum percentage 
formation was then attained. The 
percentage formation of aldehyde 
then decreased as rate of air supply 
was further increased and become 
nearly zero when a rate of air supply 
of 200 cc./min. was attained. 


2. The Effect of Mixture 
Strength at a Constant Rate 
of Air Supply of 50 cc./min. 
to the Annular Reactor 


(a) Rates of Oxygen Reacted and 
Carbon Oxides formed 
Figure 11 shows the effect of en- 
riching or weakening the mixture. A 
more pronounced negative tempera- 
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ture coefficient of reaction was ob- 
tained with the rich mixture and a 
mere complex low temperature 
mechanism was in operation for the 
weak mixture. In the latter, a maxi- 
mum was obtained during the initial 
rapid increase in oxygen consump- 
tion, after which the rate of oxygen 
consumption declined slightly, re- 
mained steady over a_ prolonged 
teperature range, and then dropped 
again until the high temperature 
reaction began. 
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Figure 8—Rate of reaction in the 

annular reactor with pentane for a 

correct mixture with air supplied at 

50 cc./min. 
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Figure 9—Rate of formation of alde- 
hydes in the annular reactor as affected 
by temperature. 
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Figures 12, 13, and Figure 8, ante 
show the variations in carbon monox- 
ide and dioxide formation with 
temperature, for various mixture 
strengths. In the rich mixture, the 
monoxide formation exceeded that of 
the dioxide over the whole range of 
temperatures as opposed to the effect 
previously described for the correct 
mixture. The weak mixture on the 
other hand, showed no formation of 
carbon oxides in the low temperature 
region. However, in the high tem- 
perature region the products were 
similar to those obtained with the 
correct mixture except that the rate 
of formation of carbon monoxide was 
less than that of the dioxide. 


(b) Production of Aldehydes and 
Peroxides 

The influence of mixture strength 
on aldehyde production at 430°C. for 

1 flow of 100 cc./min. is given by 
ri igure 14, which shows that the alde- 
hyde content of the effluent increased 
steadily with increased mixture 
strength over the range 25% weak 
to 100% rich. When calculated on 
the basis of moles of aldehyde pro- 
duced per mole of pentane however, 
a variation of from 0.16 to 0.22 only, 
was noted over this range of mixture 
strengths. 

The analyses for peroxides indi- 
cated that their rate of formation was 
negligible in the conditions of these 
and all preceding experiments with 
the annular reactor. 


Discussion 

The narrow annulus provided a 
combustion space in which the flow 
of mixture was always laminar. Thus 
reaction in the low temperature range 
began by the oxidation of the sta- 
tionary layer of mixture on the sur- 
faces and its continuance and rate 
depended upon the release of pro- 
ducts from the surface, their rate of 
diffusion into the moving stream and 
the rate of replacement of the layer 


by diffusion of unreacted mixture 
from the stream. Thus, the velocity 
of the low temperature reaction 


would be expected to depend largely 
on the time of residence in the re- 
actor, as has been demonstrated. 
The occurrence of a negative tem- 
perature coefficient of reaction over 
a subsequent temperature range can 
be attributed to the increasing ac- 
tivity of the surface in this range in 
the further reaction of an active 
intermediate product of the low 
temperature reaction. It was shown in 
Part 23 (1% that the occurrence of 
ignition in the low temperature range 
depended on the existence apart from 
the boundary layer of a central 
stream in which there were circulat- 
ing convection currents. It was in- 
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Figure 12—Rate of reaction, 50 
weak mixture supplied to the annular 
reactor at 50 cc./min. 
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Figure 14—Volume percentage of aldehyde formation as affected by mixture 
strength with annular reactor temperature maintained at 430°C. 


dicated by the experimental results 
that ignition occurred when alde- 
hydes initially formed on the surface 
were carried from the boundary layer 
into the central stream. It has been 
shown here in Figure 9 that the rate 
of formation of aldehydes rose to a 
maximum with temperature in a 
manner similar to that of the overall 
rate of oxygen consumption and then 
decreased rapidly. Since low temper- 
ature ignitions did not occur in these 
conditions it is suggested that further 
reaction of aldehydes on the surface 
at an increasing rate with temperature 
within a certain temperature range 
accounts for the negative temperature 
coefficient of the reaction. Although 
the particular aldehyde responsible 
for the effect has not been deter- 
mined, these results appear to confirm 
experiments prev iously reported in 
which the addition of acetaldehyde to 
hydrocarbon-air mixtures was shown 
to autocatalyze the first stage reac- 
tion 2). The decomposition of this 
substance to final products at in- 
creasingly higher temperatures under 
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high residence time conditions or its 
further oxidation under low residence 
time conditions accounts for the de- 
creased rate of oxygen consumption 
in this temperature range. 

The occurrence of a range of 
temperatures in which no reaction is 
apparent can be accounted for, there- 
fore, by the combination of the two 
factors given above. Thus when the 
residence time is short enough so that 
the rate of diffusion of products from 
the surface and reactants to the sur- 
face is inappreciable and the tem- 
perature is such that complete reac- 
tion of aldehydes on the _ surface 
occurs, a condition is attained such 
that no oxidation or decomposition 
products appear in the effluent and 
the overall rate of reaction is so low 
as to be nonmeasurable. This state of 
affairs continues until the tempera- 
ture is raised to a point where a new 
reaction mechanism is operative and 
where the rate of diffusion of the 
products has been raised to a value 
comparable with the residence time. 
This is not in accordance with state- 


ments by Prettre ‘!?)) and others that 
chain branching occurs on the wall 
and that deactivation of the carrier 
occurs in the gas phase and is 
accelerated by raising the tempera- 
ture. Instead, it has been shown that 
the deactivation is a surface effect 
and in consequence it is not necessary 
to devise a gas phase _ reaction 
mechanism or to make use of any of 
those devised by others in order to 
account for the negative temperature 
coefficient. 

Experiments involving more com- 
prehensive analyses of the oxidation 
and decomposition products are being 
carried out to determine the nature 
of the substances formed in the 
temperature range over which the 
coefficient of the reaction velocity is 
negative. Further experiments to 
demonstrate the powerful influence 
of surface on this phase of the reac- 
tion are also in progress. 
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Flooding Phenomena in 


Packed Gas Absorption Columns 


G. W. MINARD? and M. D. WINNING* 


This paper embodies the results of an investigation of pressure drop 
and flooding in a 6 in. diameter gas absorption column packed with 


half inch Raschig rings and Fiberglas. 


Pressure drops and flow rates 


were measured for eight combinations of gas and liquid phases in a series 


of 437 runs. 


ranges. 


The flow rates and fluid properties varied over wide 
Data were obtained on visual flooding 


This investigation attempts to show that pressure drops may be 
predieted by an empirical relationship which appears to depend on gas 


flow rate, liquid flow rate and fluid properties. 


The visual flooding 


data correlate after the method of Sherwood, Shipley and Holloway. A 
new relationship applicable to pressure drop flooding is presented in the 


Ap Dy" 
<7 bi! = 0,62 
PL a 


form 


Introduction 


Bid paper reports the results of 
an investigation which dealt with 
the relationships between flow rates 
and pressure drop, between flow 
rates and visual flooding, and between 
fluid and packing properties and 
pressure drop flooding. 

Flooding in packed absorption 
columns has been previously defined 
as: the build-up of liquid on top of 
the packing; the entrainment of liquid 
in the gas phase; the spraying of 
liquid from the top of the column; 
the rise of “slugs” of foam through 
the packing. Pressure drop flooding 
has been defined as the point, on a 
graph of pressure drop across a 
packed column versus gas flow rate, 
at which the pressure drop increases 
excessively with a slight increase in 
gas flow rate at a constant liquid flow 
rate. Both “visual” and “pressure 
drop” definitions are used as means 
of studying the data in this investiga- 


tion. Similarities and differences be- 
tween the two definitions are re- 
ported. 

In the design of packed absorption 
columns, both flooding and the pres- 
sure drop due to fluid flow through 
a packing must be considered. Flow 
rates usually set the column diameter. 
Where these flow rates have been 
established by an economic study, 
they must be checked to insure that 
they will not produce flooding with 
its resultant loss of liquid and build- 
up of back pressure on the gas stream. 
When blowers are being sized the 
pressure drop across the column 
must be available. Flooding can 
seldom be observed visually, but the 
pressure drop across a column may 
be observed on a differential pressure 
gauge. Thus the prediction of the 
pressure drop below which flooding 
is unlikely to occur is useful to 
operating personnel. 


Both visual flooding and pressure 
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drop are related to flow rates. For 
visual flooding, liquid and gas ve- 
locity relationships have been given 
by Sherwood, Shipley and Hollo- 
way“), The relationship between 
pressure drop and flow rate is usually 
depicted graphically. Early  inter- 
pretations ‘?- 3) of the resulting logar- 
ithmic graphs postulated two inter- 
secting straight lines with a definite 
break point or discontinuity as shown 
in Figure 1. Attempts were made by 
various investigators to relate this 
discontinuity to visual flooding. 
Zenz ‘), on the other hand, sug- 
gested that the graph should not show 
a discontinuity, but rather continuous 
pressure build-up with increasing 
flow rate. Data obtained by Ufford °°) 
bore out this contention, and the 
authors on examining the data in the 
literature found that the two straight 
lines could be replaced by a smooth 
curve. This paper suggests a form of 
empirical equation to fit the proposed 
smooth curves. 

In the pressure drop flooding de- 
finition the flooding point would be 
that at which the tangent to the 
pressure drop-flow rate curve is a 
vertical line parallel to the pressure 
drop axis. Since much of the data in 
the literature does not reach this 
region, it was found convenient to 
take the pressure drop at flooding as 
that pressure drop at which the slope 
of the tangent to the pressure drop- 
flow rate curve exceeded 85° 

The authors conclude from the 
data examined (3.4.5) that this par- 
ticular pressure drop associated with 
flooding in the pressure drop flooding 
definition, appears to be independent 
of flow rates, and is instead a func- 
tion of packing and fluid properties. 
The points marked (*) in Figure 1 
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° 
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Figure 1—Typical pressure drop-flow rate curve. 


are considered to fit the definition 
of pressure drop flooding. It may be 
seen that the pressure drop at these 
points is constant, being independent 
of the liquid and gas flow rate. The 
nature of the relationship between 
pressure drop and packing and fluid 
properties will be discussed in this 


paper, Pr, 
Equipment 
Basically, the equipment consisted GAS HOLDER GAS 


FLOWMETERS 
of a packed column common to two 


flow circuits. The liquid flow circuit 
consisted of the column, a centrifugal 
pump, a supply tank, a constant head 
tank, a pair of rotameters, and the 


pipe, fittings, and valves necessary for 
connecting the parts and for control- 


ling the flow rate. The gas flow cir- 
cuit consisted of the column, a sLoweR| 
blower, a gasholder, a pair of rota- 
meters and the necessary pipe, fittings, 


and valves. 


° . . TOR 
The instrumentation was relatively Ls 0 


simple. Two double clean-out mano- 
meters were connected to the column 
via a valve-manifold system to per- 
mit measurement of the differential 
pressure across the packing. Static 
pressures were also measured on these 
manometers at the gas meters, gas- 
holder, column base and column top. 
Copper-constantan thermocouples 
were used, in conjunction with a 


Symbols V—valve 


potentiometer, to measure tempera- 
tures, 


Figure 2 is a schematic diagram of 
the equipment. 


The Column and Packings 

The column was a standard 6 fet. 
length of 6 in. I.D. Pyrex glass pipe 
mounted vertically in a supporting 
framework. The ends were attached 
to metal plates by aluminum flanges. 
Further connections were then made 
to the metal plates which also served 
as a means of attaching the pipe to 
the rigid metal framework. The pack- 
ing support was attached to the 
underside of the lower plate. This 
support consisted of one-quarter in. 
wire mesh basket 4 in. deep and 8 in. 
square. This type of support was 
chosen to prevent any reduction of 
flow area at the support. 

Two column packings were used. 
One-half in. Raschig rings were em- 
ploy ed in the first part of the investi- 
gation, and “Fiberglas” was employed 
in the second part. The rings were 
dumped into the water-filled column 
to give a packed height of 5 ft. The 
fiberous packing was made by cutting 
6 in. circles from one inch “Fiberglas 
Dustop Air Filters”. The short cylin- 
ders so obtained were then pushed 


HEAD 
TANK 







LIQUID 
FLOWMETERS 


Tececeey i 


SUPPLY 
av TANK 
COLUMN 


Figure 2—Equipment flow sheet. 


T—thermocouple well 
Pct—pressure tap at column top 
Pcb—pressure tap at column bottom 
Pri—pressure tap at large gas rotameter 
Prs—pressure tap at small gas rotameter 
Ph—pressure tap at gasholder. 
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up inside the dry column from the 
bottom. A variation in the fraction of 
frec void was obtained by varying 
the height of a constant mass of 
fiberous packing. 

A seven-holed liquid distributor at 
the top of the column provided uni- 


forin liquid distribution. at 
N 

Column Operation for Flooding 2 
and Hold-up Measurements 2 
To determine flooding conditions a 
in the packed column, gas and liquid = 
were circulated at constant rates. 3 
When stable operation had been ob- ot 
tained at one particular set of flow y 
rates, readings were taken of the 3 
pressure drop across the column, of ” 
gas and liquid temperatures, of the x 


pressure in the gas rotameter and of 
room temperature at the panel board. 
Periodically, during the runs, the 
barometric pressure was noted. 


After recording the flooding data 
the dynamic liquid hold-up was de- 
termined in the following manner. 
While the column was operating, the 
height of liquid in the bottom section 
was noted. Then all flows including 


t 
the column discharge were simul- 





ea O. 
tancously shut off. After a draining o 100 1000 
period of ten minutes the liquid level GAS FLOW RATE "G" LOS 4, ey2 
Figure 4—Pressure drop vs. flow rate for the system carbon dioxide-water flowing 
LBS 
through Raschig rings. The parameter is the liquid flow rate “L” HR. FT? 


in the bottom section was again re- 
corded. The difference of the two 
readings was taken as the dynamic 
liquid hold-up. 

The above procedure was then 
duplicated for other flow rates. 

For those systems other than air- 
water- Raschig rings no hold-up data 
were collected. Several runs were 
made with the packing wet but with 
no liquid flowing. At the conclusion 
of those runs with wet packing, the 
top section of the column was re- 
moved, and the packing was dried by 
blowing air through it. Then pres- 
sure drop- -flow rate data for the dry 
packing were secured. 


Hof T_ 


Experimental Results 





Visual Observations 
1. Column Operation with Raschig 
Ring Packing 


PRESSURE OROP AP INS. 


° 


The flow characteristics observed 
for Raschig ring packing agreed 
closely with those of White (*) 

At low liquid rates the flow of 
liquid through the packed column 
appeared to vary with gas flow rate. 
At low gas rates the liquid formed 

drops on the packing, and as these 

10 100 +9" os 0009 became larger they flowed dow 
GAS FLOW RATE "G" LBS/.p p72 vecame larger they flowe own 

Figure 3—Pressure drop-flow rate curves for the system air or nitrogen and from one unit of packing to the next 


water flowing through Raschig rings. The parameter is liquid flow rate below, For a higher gas rate, the 
“L” LBS/HR. FT?; air ©; nitrogen @. liquid hold-up in the column had 
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Figure 5—Pressure drop-flow rate curves for the system air-carbon tetrachloride 


flowing through Raschig rings ,»L = 99.5 LBS/FT*. The parameter is liquid 
flow rate “L” LBS/HR.FT? 
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Figure 6—Pressure drop-flow rate curves for the system air-carbon tetrachloride 

and ethanol flowing through Raschig rings pL = 74.6 LBS/FT*. The parameter 
is liquid flow rate “L” LBS/HR.FT2 


accumulated enough to form con- 
tinuous streams instead of drops, 
Further increase in gas flow rate in- 
creased the size of the streams until 
the liquid appeared to become the 
continuous phase with the gas bub- 
bling through it. Generally at this set 
of conditions a layer of liquid or 
foam formed on the top of the pack- 
ing. 

At high liquid rates (say L ~ 
10,000 Ibs./hr./ft.2) the first liquid 
flow was not in the form of drops but 
rather in continuous streams; i.e., the 
first flow stage mentioned in the 
paragraph above was not present. In- 
creasing the gas flow rate had the 
same effect as previously noted. 
the higher liquid rates the layer of 
foam noted at flooding became a 
layer of liquid with only a few gas 
bubbles in it. 


The behaviour of the column fol- 
lowed the same general pattern for 
variations in gas and liquid rate. 
When the flow rates are high enough 
to cause visual flooding, a slight dif- 
ference is noted in the lay er of liquid 
on top of the packing. 


2. Column Operation with 

Fiberglas Packing 

The effect of gas and liquid flow 
rates on the appearance of the Fiber- 
glas packed column did not differ 
greatly from that for ring packing 
provided the flow rates were below 
those which would cause flooding. A 
few major differences in appearance 
were noted between the ring packed 
and Fiberglas packed columns when 
visual flooding had been achieved. 


As flooding conditions were ap- 
proached the packing appeared to fill 
completely with liquid. However, a 
slight increase in gas velocity beyond 
flooding caused the packing to loosen 
or expand, and to move about in the 
column. The expansion opened up a 
number of large voids in the packing. 
A further slight increase in gas ve- 
locity caused a build-up of liquid on 
the top of the packing as well as an 
increase in the size of the void. 


A greater increase in gas velocity 
over that for flooding caused such a 
loosening effect that channelling de- 
veloped, and the liquid head on top 
of the packing disappeared. An ad- 
ditional increase in gas velocity only 
aggravated this condition. When the 
gas and liquid flows were stopped, 
after a period of flooding, it was 
found that the top layer of fiber had 
been broken into short lengths by 
the agitation of the liquid. Under 
further operation some of these fibers 
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could be noted being blown up to the 
top of the column. 

‘the loosening and breaking of the 
Fiberglas packing resulted in a change 
in the fraction of free void. This 
change could not be measured. 

‘The operation of the column was 
stable for all flow rates below, at and 
above those required for visual flood- 
ing in both packings used. For opera- 
tion at flooding a period of time was 
in the 


necessary for the conditions 
column to become steady. However, 
a stable operation was always 


achieved. 


Correlation of Results 


1. Pressure-Drop and Flow Rate 

Relationships 

a. Raschig Ring Packing: The first 
data collected were obtained for the 
purpose of finding the variation in 
pressure drop-flow rate relationships 
due to variations in gas properties. 

Water was used as the liquid and 
air, nitrogen, and carbon dioxide 
were used respectively as the gas. 
The data for air and nitrogen fell so 
closely together that the same graph 
could be used to represent both. The 
data for carbon dioxide did show 
significant differences. Changes in gas 
viscosity are rather small so that the 
effects shown were generally due to 
gas density. 

Figure 3 presents the relationship 
between pressure drop (AP) and 
gas flow rate (G) with liquid flow 
rate (L) as the parameter, for the 
systems air-water-Raschig rings, and 
nitrogen-water-Raschig rings. Also 
included are the base lines for air 
flowing through dry packing and 
moistened packing. Figure 4 presents 
the corresponding data for the sys- 
tem carbon dioxide-water-Raschig 
rings. 

The data from which the Figures 
3 to 11 have been prepared are tabu- 
lated in the Ph.D. thesis ‘®), by M. D. 
Winning entitled, “Flooding Pheno- 
mena in Gas Absorption Columns”. 
1952, University of Toronto, which 
is available from the University 
library. ’ 

Several important features are im- 
mediately evident from Figures 3 and 
4. First it will be noted that the 
data fall on smooth curves. This bears 
out the contention of Zenz ), It is 
true that for the lowest liquid flow 
rate, a portion of the curve appears 
straight, but the curvature is very 
evident at the higher liquid rates. it 
will also be noted that the curves 
straighten out at the higher liquid 
rates. 


Although it is difficult to distin- 
guish between the data for air and 
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Figure 7—Pressure drop-flow rate curves for the system air-carbon tetrachloride 
and ethanol flowing through Raschig rings pL = 61.5 LBS/FT*. The parameter 
is liquid flow rate ““L” LBS/HR.FT? 
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Figure 8—Pressure drop-flow rate curves for the system air-carbon tetrachloride 
and ethanol flowing through Raschig rings pL = 56.8 LBS/FT*. The parameter 
is liquid flow rate “L” LBS/HR.FT? 
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nitrogen, comparison of these with 
the lines for carbon dioxide show 
that increasing the gas density in- 
creases the amount of the gas passed 
at the same pressure drop. 


As a result of the fact that the 
curves appear to be continuous, at- 
tempts were made to find empirical 
equations to fit the curves. Equations 
for gas flow through dry or wet 
packings were easily obtainable in 
the form AP = q Gr where AP is 
the pressure drop in inches of water 
per foot of packing height, “G” is 
the gas flow rate in pounds per Be 
square feet of column area, and “q” 
and “r’’ are constants. When gas and 
liquid are present in the column in 
counterflow the pressure drop ap- 
peared to increase more rapidly than 
would be indicated by a simple 
power function of the gas flow rate. 


Of a number of equations tested 
against the data the relationship 


sufficient data in the literature, it has 
not been possible to obtain a general- 
ized correlation between them. 


A second set of data was obtained 
for mixtures of ethanol and carbon 
tetrachloride as liquid, air as gas and 
for Raschig ring packing. Mixtures 
of these two liquids allow a wide 
variation in liquid density with but 
slight change in viscosity. Figures 5 
to 9 show the effect of variation of 
liquid density on the pressure drop- 
flow rate relation. 


It may be noted that changes in 
liquid density have two effects. First 
at one liquid velocity, the gas 
capacity of the column increases, with 
increasing liquid density, as shown 
by a shift of the curves to the right. 
Second, the liquid capacity of the 
column is increased as shown by the 
great variation in the positions of the 
curves at high liquid flow rates. Just 
what fraction of the gas capacity 
increase is due to liquid density 


1G changes and what part is due to gas 
AP = qG' + stan{ - : density changes cannot be deter- 
a inined. 


was found to fit the data reasonably 
well. The terms q, r, s, and t were 
found to vary with liquid flow rate 
and gas density, but because of in- 


As with the first set of data, at- 
tempts were made to obtain empirical 
equations. In all cases it was found 
possible to relate the pressure drop 


TABLE 1 
: Ap p\* 
CaLcuLatep VaLuEs or — by! 
prDe \a 

i. Schoenborn........0...5 ..0.630 8. Elgin and Weiss“ ......... 0.542 
De kev edvvewenesas ies : nee ae 9. Elgin and Weiss...........0.578 
Bee athictie wwsia's (op ocs wate hci 10. Elgin and Weiss...........0.648 
4. NE ace hn ea nla es Us RNR! oo ee ci io ke os 0.648 
ares ee er RMS (koe cA KGa: Caos OR 0.605 
De ie Shee Cp iue 0i .. .0.602 1S.. Awthore......... . .0.667 
7. Dougherty®...............0.645 14. Authors...... .0.622 

TABLE 2 


PackING & Gas/Liguip SysTEMS USED IN PRESSURE Drop FLOODING CORRELATION 


Investigator Pac king Type & Size Gas Liquid 
Schoenborn and Dougherty. ‘ Carbon Rings ” Air Water 
C eramic Rings ” 
Porcelain Rings 14” | and 
Berl Saddles 14” 
Carbon Rings yr | Oil 
Elgin and Weiss“... . « .| Berl Saddles 14” Air Water 
| Berl Saddles 14” 
| Clay Rings 5B” | 
| Clay Balls 4” 
| 
BOs: 2 cies Sie Rew clerra es Rings 114” Air Water 
| Rings 1” 
| Rings 34” 
Rings 34” 
| Rings 14” 
Rings +” 
Peis ce een Reviewed all the above plus 
Ceramic Rings i” | Air Water 


and the gas flow rate by the equation 


G 
AP = qG' +s tan ( ) 
2 


The influence of free void fraction 
was not investigated nor was the 
influence of packing orientation. 

b. Fiberglas Packing: Fiberglas has 
been suggested as a possible packing 
material for columns on the basis of 
the large surface area obtainable per 
cubic foot of packed volume. In 
order to improve on the other pack- 
ings already in use, an arrangement 
of packing to give an area greater 
than 250 sq. ft per cu. ft. is suggested. 
Tests were made on two packings of 
Fiberglas having a calculated area of 
270 and 412 sq. ft. per cu. ft. respec- 
tively. The fractions of free void for 
these packings when wet were 0.95 
and 0.91 respectively. 

The tests on the packing having a 
free void fraction of 0.91 indicated 
that the packing would pass gas 
without difficulty, but flooding oc- 
curred even at the lowest gas rates 
when there was any liquid flow. The 
runs made with the packing having a 
free void fraction of 0.95 indicated 
allowable flow rates high enough to 
warrant using the packing indus- 
trially. The data obtained are shown 
in Figure 10. 


The data for Fiberglas packing 
show that much smaller pressure 
drops occur than for other packings 
with equivalent flow rates. Thus, a 
tower packed with Fiberglas would 
have a greater capacity, both for 
liquid and gas, than a tower filled 


with the more common. industrial 
packings. However, when the _in- 
crease in pressure drop which 


characterizes flooding, occurs, it is 
extremely abrupt. The lifting or open- 
ing of the packing mentioned under 
visual observations creates some doubt 
as to whether the pressure drop re- 
corded is the result of flooding or due 
to lifting of the packing. 


The doubt as to the reliability of 
the fraction of free void at flooding, 
together with the knowledge that the 
Fiberglas disintegrates under the agi- 
tation of flooding, restricts the range 
of data for which empirical equations 
may be determined. However, for the 
flow of air and water through 
Fiberglas having a fraction of free 
void of 0.495, where the gas rate does 
not exceed 700 lbs./hr./ft.2 and the 
liquid rate does not exceed 10,000 
lbs./hr./ft.*, pressure drop (AP) may 
be evaluated as 


AP = (39,500 + 2.42L) x 10-9 G4 


L = Liquid flow rate in pounds per 
hour square feet of column area. 
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‘Yhe restrictions applied here still 
allew a range of operating variables 
which should be of use industrially. 
However, no attempts should be 
made to use Fiberglas near flooding 
conditions without introducing some 
means of restricting its movement. 


2. Flooding Correlations 


Visual Flooding: The correlation 
method of Sherwood, Shipley and 
Holloway “ for visual flooding flow 
rates, was applied to the data of this 
investigation as a means of comparing 
them with those of other investiga- 
tors. The runs which fit the definition 
of visual flooding were plotted as 


L /pg  Ugtapg ,, 
= Vs. a 
G PL geé® pi 


and are shown in Figure 11. The solid 
line represents the correlation pro- 
posed by Sherwood et al., as the best 
locus for all previous data. The two 
dashed lines indicate the width of the 
band covered by previous data. It 
can be seen that most of the runs of 
the present investigation fit the pro- 
posed correlation remarkably well, 
and all fall within the indicated band 
except one of the doubtful runs using 
Fiberglas. 

This shows that the construction of 
equipment, collection of data and 
methods of calculation are com- 
parable to those of previous investi- 
gations and that the comparison of 
data with those previously reported 
is justifiable. 

Pressure Drop Flooding: On ex- 
amining the data of the flow region 
where “the pressure drop appears to 
increase without limit, certain fac- 
tors seemed to be related. 

These factors were: 

Pressure drop—Ap_pounds/square 
foot 

Packing size (diameter)—D feet, out- 
side diameter for spheres, fibers, 
and rings (usually equal ring 
length), face to face distance for 
saddles 

Packing void fraction—¢ cu. ft. of 
free space/cu. ft. of packed volume 

Packed surface area—a sq. ft. of 
area/cu. ft. of packing 

Liquid viscosity—u,, pounds/foot sec. 

Liquid density—P , pounds/cu. ft. 

Dimensional “a sis was applied to 
these factors, and the followi ing 
equation resulted: 


dp _[(a) (Dixy | 
Dpt D 5 Mu? = 


Various arrangements of these 
dimensionless ratios were tried as a 
means of correlating the pressure 
drop at flooding. Values of the pack- 
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Figure 9—Pressure drop-flow rate curves for the system air-carbon tetrachloride 


and ethanol flowing through Raschig rings pL = 53.4 LBS/FT°. The parameter 
is liquid flow rate “L”? LBS/HR.FT? 
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Figure 10—Pressure drop-flow rate curves for the system air and water flowing 
through Fiberglas. The parameter is liquid flow rate “L” LBS/HR.FT?; 
e« = 0.950; ¢ = 0.9le. 
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Figure 11—Visual flooding correlation after Sherwood. Raschig rings O Fiber- 


glas @. Solid line represents Sherwood’s 


correlation. Dashed lines indicate range 


of data used in Sherwood’s correlation. 


ing and liquid properties, 
with the pressure drop were secured 
from the literature for 14 different 
sets of conditions. In four of these, the 
liquid properties were not reported, 


together 


and it was necessary to assume 
values. The assumed values were 
taken as the appropriate liquid pro- 
perty at 70°F. and one atmosphere 
pressure. As a result, the possible 
error in the calculated dimensionless 
ratios is quite large. 


When values of 


Ap 1/3 
2 ’) 


were calculated, they were found to 
be almost constant. However, a 
variation with liquid viscosity was 
noticed. Finally, multiply: ing by pty, 
where uy, is in centipoise instead of 
pounds per foot second, gave 


Ap /dD\ 
piDe \a oa 


Table 1 gives values 


Ap D ee 0.1 
prDe \a - 


obtained from the published data of 

other workers along with values from 

the data obtained in this investigation. 

Table 2 lists the packings and pack- 

ing size used by the other workers. 
The fact that 


Ap D PS 0.1 ? 
pl De a a" = 0.62 


means that it is possible to predict 
the pressure drop which will occur 
at the onset of pressure-drop flood- 
ing for any packing and fluid com- 
bination. Thus, for design purposes, 
it is possible to calculate the pressure 
drop which might occur in a packed 
column operating at any percentage 
of flooding capacity. 


Conclusions 


In this investigation an attempt has 
been made to study the relationship 
between pressure drop, gas and liquid 
flow rates and flooding in a gas ab- 
sorption column. It has been found 
possible to fit an empirical equation 
to flooding data, and it is suggested 
that additional data be secured in an 
attempt to relate the terms of this 
equation to packing and fluid pro- 
perties. A method of predicting pres- 
sure drop at the onset of flooding is 
suggested. 


Nomenclature 


a packing surface area (ft.?/ft.*) 

D packing diameter (ft. ) 

G gas flow rate (Ibs./hr. ft.) 

g gravitational constant (ft./sec.?) 
liquid flow rate (Ibs./hr. ft.) 

Ap pressure drop (Ibs. 
in.?) 


ft.2 or Ibs./ 


arbitrary constants 

gas velocity (ft./sec 

packing void 

liquid density (Ibs./ft.*) 

liquid viscosity (Ibs./ft. sec 
pressure drop (ins. HO : of 
packing) 
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